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Abstract 



Deep-inelastic ep scattering data, taken with the HI detector at HERA, are used to 
study event shape variables over a large range of "relevant energy" Q between 7GeV 
and 100 GeV. Previously published analyses on thrust, jet broadening, jet mass and 
C parameter are substantially refined and updated; differential two-jet rates treated as 
event shapes are presented for the first time. 

The Q dependence of the mean values is fit to second order calculations of perturbative 
QCD applying power law corrections proportional to 1/Q P to account for hadronization 
effects. The concept of these power corrections is tested by a systematic investigation in 
terms of a non-perturbative parameter and the strong coupling constant. 



Kurzfassung 

Ereignisformvariablen in der tiefinelastischen ep-Streuung, gewonnen aus Daten des Hl- 
Detektors bei HERA, werden fiber einen grofien Bereich der „relevanten Energieskala" 
Q von 7 GeV bis 100 GeV untersucht. Bereits publizierte Studien zu thrust, Jetbreite, 
Jetmasse und C Parameter werden erheblich verbessert und aktualisiert; Ergebnisse zu 
differentiellen Zweijetraten als Ereignisformen werden erstmals vorgestellt. 

Bei der Anpassung der Q-Abhangigkeit der Mittelwerte an Berechnungen der per- 
turbativen QCD in zweiter Ordnung sollen potenzartige Korrekturterme der Form 1/Q P 
Hadronisierungseffekte berucksichtigen. Eine Uberprufung des Konzepts solcher Poten- 
zkorrekturen erfolgt im Rahmen der Bestimmung eines nichtperturbativen Parameters 
57p_ 1 und der Kopplungskonstanten der starken Wechselwirkung. 
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Chapter 1 
Introduction 



Curiosity is one of mankind's elementary driving forces — the more so when scientists 
are concerned. In European history the appearance of the first "professional scientists" 
is usually dated back to antique Greece, where, besides others, philosophers laid the 
foundations for a great part of European culture. Two of them living around 500 B.C., 
Leukippos and Demokritos, suggested first that matter consists of tiny indivisible particles 
called atoms. They may be considered the forefathers of today's elementary particle 
physicists. 

Their reasoning was of a basically theoretical nature. Turning to experiment, a "quan- 
tum leap" leads us to Antoni van Leeuwenhoek (1632-1723). Employing microscopes he 
constructed himself, he was able to observe "creatures that swarm and multiply in a drop 
of water"[] for the first time. 

With the advent of the modern microscope Hadron-Elektron-Ring-Anlage HERA 
at the Deutsches Elektronen-SYnchrotron DESY in Hamburg three centuries later, the 
resolution power could be increased by about twelve orders of magnitude. Now it is 
possible to study "gluons that propagate and split in a droplet of nuclear matter called 
proton." 

The theory that is believed to describe the strong interaction, responsible for the struc- 
ture of hadronic (nuclear) matter, is called Quantum ChromoDynamics (QCD) in analogy 
to Quantum ElectroDynamics (QED) dealing with the electromagnetic interaction. Both 
are part of a more comprehensive theory usually referred to as the Standard Model (for 
introductory textbooks s. e.g. |2|, |3|, 11). It unites three, the weak, electromagnetic and 
strong interaction, out of the four fundamental forces observed in nature into a common 
framework. Gravity so far does not fit in. Comprising merely twelve elementary particles 
of spin 1/2, six quarks and six leptons, and the spin 1 exchange quanta of the three forces, 
i.e. the Z° and W ± bosons, the photon 7 and eight gluons g, the Standard Model gives a 
very good account of a tremendous amount of data gathered by experiments during the 
last decades ||. 

The HERA collider offers the unique possibility to investigate the structure of the 
proton and its interacting constituents in a completely new kinematic domain unreached 

1 H.G. Wells, The War of the Worlds, 1898. 
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by fixed-target experiments. With respect to this analysis, the large accessible range in 
"relevant energy" Q is of great advantage. Observables can that way be studied and 
compared to theory in dependence of the available energy Q in a single experiment. To 
that goal, the thesis is organized as follows: 

Chapter two describes the machinery necessary to perform ep scattering experiments 
at the required high energies with special attention to the parts relevant for this analysis. 
The third chapter explains the basic kinematics of the scattering process. Subsequently, 
the observables to measure are defined in chapter four. The series of chapters five, six and 
seven deals in detail with the selection of data, their comparison with simulations and 
the correction for detector imperfections. Fully corrected data are presented. The two 
following chapters eight and nine are dedicated to the derivation of the necessary theoret- 
ical input employing perturbative QCD (pQCD). Finally, experimental and theoretical 
results are compared and combined to test the power correction approach to hadroniza- 
tion effects that was initiated by Dokshitzer and Webber ||. The last chapter presents a 
summary and an outlook. 



Chapter 2 
HERA and HI 



2.1 The Storage Ring HERA 

Scattering experiments are the basic tools to investigate the elementary building blocks 
of matter, i.e. quarks and leptons, and the fundamental forces acting between them. 
Two rather straightforward set-ups are e + e~ and pp colliders. In the last decades, a lot 
of valuable information has been accumulated with such machines |J. However, for the 
detailed study of a strongly bound, hadronic object like the proton, it is preferable to use a 
probe that does not itself interact strongly. Fixed-target experiments employ electrons and 
myons as well as neutrinos to determine the structure of nucleons. To enlarge the resolving 
power, it is necessary to increase the energy available in the collisions by accelerating the 
"targets," i.e. a collider experiment is asked for. 

The Hadron-Elektron-Ring-Anlage HERA, shown in fig. |2.1| , is situated at DESY 
in Hamburg, Germany, and is currently the world's only facility where electrons^ and 
protons are accelerated in two separate storage rings to final energies of 27.5 GeV and 
820 GeV respectively. The resulting center of mass energy of about 300 GeV corresponds 
to electron beams of 50 TeV for fixed-target experiments. For a recent overview of the 
knowledge gained on the structure of nucleons consult e.g. |7J. 

Along the circumference of 6.3 km two locations, the north hall and the south hall, are 
assigned to the study of ep scattering. The interaction regions, where during operation 
every 96 ns particle bunches may collide, are surrounded almost hermetically by complex 
detectors, HI || and ZEUS . They are dedicated to the task of measuring as many col- 
lision products as precisely as possible. Two other experiments, HERMES and HERA-B, 
are situated in the east and west halls and are committed to spin and B meson physics 
respectively. 



1 Due to a considerably higher lifetime at large currents, the electron beam has been replaced by a 
positron beam in July 1994. Since for the purpose of this study it does not make a difference, the term 
electron will henceforth be used synonymously for positrons, too. 
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Figure 2.1: The storage ring HERA at the DESY laboratory in Hamburg, Germany. 

2.2 The HI Detector 

The basis of this experimental analysis are data collected with the HI detector, which is 
located at the northern interaction point of HERA. Only a brief overview of the total 



system, shown in fig. |2.2| , will be given in the next section. The parts that provide the 
main information needed in this study will subsequently be discussed in more detail. 
A complete description of the detector and its performance in the first three years of 
operation (1992-1994) can be found in [g gg. 

During the winter shutdown 1994/95, a major upgrade of the HI backward region 



was undertaken [IT]]. However, none of the components relevant to this analysis were 
significantly affected. Therefore, these improvements will not be covered here, although 
data from 1994 up to 1997 are used. Events where the scattered electron was found in 
the backward calorimeter are taken from 1994 data only. 



2.2.1 Overview of the HI Detector 



The most obvious feature of the HI detector in fig. |2.2| is its asymmetric design. Because 
the momentum of the protons is about 30 times higher than that of the electrons, the ep 
center of mass system is moving along the direction of the proton beam. As a result, the 
density of collision products hitting the equipment in the "forward" region is very high 
which is accordingly reflected in the more massive as well as finer granulated material 
of that part. The term "forward" refers to the conventional coordinate system used in 
the HI collaboration where the proton beam direction is defined to be the +z-axis. The 
x-axis (y-axis) points from the nominal interaction point at the center of the HERA 
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ring (upwards). Subsequently, the detector components are briefly described proceeding 
from the innermost parts outwards: 

• The tracking system: 

Directly surrounding the beam pipe and beam magnets, a tracking system is in- 
stalled to measure the momentum of charged particles. It consists of two main 
parts: the Central Tracking Device (CTD) covering the region around the nominal 
interaction vertex and the Forward Tracking Device (FTD) supplementing it in the 
+2;-direction. 

• The calorimeters: 

To complement the momentum measurement and to detect neutral particles, the 
CTD and FTD are enclosed in the forward and central region by a sampling calo- 
rimeter (LAr) with Liquid Argon as active material. For the innermost absorber 
stacks, lead has been chosen to ensure a good containment and energy determina- 
tion of electromagnetic showers produced by electrons and photons (Electromagnetic 
CALorimeter, ECAL). The outer part of the LAr (Hadronic CALorimeter, HCAL) 
predominantly measures hadronic showers and is equipped with steel absorber 
plates, which also serve as mechanical support structure. 

The remaining holes of the LAr around the beam pipe are closed with a silicon- 
copper calorimeter for polar angles below 4° (PLUG) and a lead-scintillator 
calorimeter in the backward direction (Backward ElectroMagnetic Calorimeter, 
BEMC). During the upgrade in 1994/95 the BEMC was replaced by a lead- 
scintillating fibre Spaghetti Calorimeter (SpaCal) ||12|| , which is subdivided into 
a first electromagnetic section and an additional second part to determine energy 
leakage and to improve the containment of hadronic showers. For the purpose of 
this analysis, it was only employed to supplement eventual energy deposits in the 
backward direction not due to the scattered electron. 

• The superconducting coil: 

A cylindrical superconducting coil providing the magnetic field of 1.15 T for the 
trackers envelops the calorimeters. Thereby, the amount of dead material in front 
of the calorimeters is reduced and the time of flight of myons within the magnetic 
field is increased improving the resolution of their momentum measurement. 

• The instrumented iron yoke and myon system: 

The IRON return yoke (IRON) of the magnet, enclosing almost completely all 
other parts of the detector, is sandwiched with streamer tubes for the measure- 
ment of myons and energy leakage from the inner calorimeters (Tail Catcher, TC). 
Supplementary chambers inside the IRON further improve the evaluation of myon 
tracks. Myons with high momenta in the forward direction are analyzed by a spec- 
trometer consisting of four drift chambers in the magnetic field of 1.5 T of a toroidal 
coil in front of the IRON. 
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Figure 2.2: Schematic layout of the HI detector. 
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Figure 2.3: The layout of the HI luminosity system 



The luminosity system: 



The luminosity system sketched in fig. |2]^ utilizes the Bethe-Heitler process ep — > 
ep7- The electrons and photons are detected in coincidence by a hodoscope of 
crystal Cherenkov counters. Whereas the electrons from these grazing collisions, 
i.e. 9 e > ~ 180°, are deflected by magnets from the beam line to hit the Electron 
Tagger (ET) at z = —33.4 m, this is not possible for the photons, of course. They 
leave the proton beam through a window at z = —92.3 m instead, where the beam 
pipe bends upwards, and reach the Photon Detector (PD) at z — —102.9 m. 



The main background is caused by bremsstrahlung from residual gas atoms eA — > 
eAy. The technique of electron pilot bunches that do not have colliding proton 
counterparts allows to correct for it. 
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2.2.2 Details on the Main Components 

The basic ingredient for this analysis is the calorimetric information provided by the LAr 
and BEMC in the form of electromagnetic and hadronic clusters reconstructed from the 
primary energy depositions. The fragmentation region of the proton remnant, which is 
related to "soft" physics, can be excluded more easily after performing a Lorentz trans- 
formation into the Breit frame (s. section |4.1| ). This procedure requires a good electron 
identification. On the one hand, the electron cluster has to be separated from the ha- 
dronic final state, and on the other hand, it is necessary for the extraction of the event 
kinematics which is crucial in the calculation of the boost to the Breit frame. There- 
fore, additional information of the tracking system, especially the CTD, is exploited to 
determine the properties of the scattered electron and the actual event vertex. 

The main attention in this analysis rests on mean values and normalized distributions. 
Hence, the luminosity system needed to measure absolute cross sections is not of prevalent 
importance. 

Concluding, the CTD, LAr and BEMC are the most important subdetectors for the 
observables considered and will therefore be described in more detail JTOf: 



The central tracking device CTD: 

The CTD is subdivided into two major units, the Central Jet Chambers CJC1 



and CJC2 as shown in fig. |27| in a view perpendicular to the beam. The track 
reconstruction in the central region, 25° < 9 < 155°, depends primarily on these 
concentric drift chambers with wires strung parallel to the beam pipe. They provide 
a good space point resolution in the (r, 0)-plane of 170 /zm and, by comparing the 
signals at both ends, furthermore supply information about the ^-coordinate of the 
hits. To avoid completely insensible paths for almost straight tracks, the sequence 
of the sense wires is inclined with respect to the radial direction. 

The CJC1 is sandwiched between the Central Inner and Outer Z-chambers CIZ and 
COZ with wires oriented in the (r, 0)-plane. They complement the measurement 
in the CJC's with data on the z-coordinate of track elements with a precision of 
^ 300 (Ltm. 

Mainly for the purpose of fast timing (< 96 ns) and triggering, multiwire propor- 
tional chambers are added for polar angles between 5° and 175°. Two of them, the 
Central Inner and Outer Proportional chambers CIP and COP, are indicated in 



fig. [2.4| . For electrons in the backward direction, only short track segments can be 
seen in the CTD. The Backward Proportional Chamber BPC, consisting of four 
sensible layers of wires oriented along azimuthal angles of 0°, 45°, 90° and 135°, sub- 
stantially improves the tracking for polar angles of 155.5° < 9 < 174.5°. Installed 
directly in front of the BEMC, fig. |2.5| , it also helps to distinguish between photons 
and electrons in this region. 

In general, the tracking system was designed to determine the momentum and angles 
of charged particles to a precision of a p /p 2 ~ 0.003 GeV -1 and er e ~ 1 mrad. 
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The liquid argon calorimeter LAr: 

As mentioned, the LAr, covering polar angles of 4° < 9 < 154°, is a sampling 
calorimeter where layers of absorber material and sensitive gaps, filled with liquid 
argon, alternate. Since only a sample of all energy deposits can be measured that 
way, the energy resolution of this technique is worse than that of calorimeters built 
of only one medium. To achieve nevertheless a maximal precision, care has been 



taken that the orientation of the absorber plates, shown in the upper part of fig. |2.6| , 
ensures angles near the normal direction of the stacks for particles originating from 
the vertex. Another drawback is the difference in response on electromagnetic and 
hadronic showers produced by impacts of particles of the same original energy. The 
LAr is n on- compensating. 

This is counterbalanced by stable calibration and homogeneous response properties 
as well as the possibility of a compact but still finely segmented construction. The 
structure of the read-out cells, again in an (r, z)-view, is presented in the lower half 
of fig. 

Apart from the general subdivision in an Electromagnetic and a Hadronic section, 
the LAr is composed of eight wheels each further partitioned into eight octants 
in (f). Due to the use of weak lead for the ECAL, a pointing geometry, blind for 
hits straight on the border between two octants, could not be avoided in contrast 
to the CTD. Particles reaching these ^-cracks are detected only in the non-pointing 
HCAL. All wheels are, with respect to the nominal interaction point, combined to 
four aggregates: a Backward, Central and Forward Barrel and the Forward end 



cap with an Inner and Outer ring. The abbreviations in fig. ^T6| translate e.g. as 
"outer hadronic ring of forward end cap wheel 1" for OF1H. 

Test beam results yield energy resolutions of ~ 11%/ ' \[E © 1% for electromagnetic 
and ~ 50%/\/E © 2% for hadronic showers.^ The absolute energy scales are known 
to 1-3% |13j and 4% for electrons and hadrons respectively. 



The backward electromagnetic calorimeter BEMC: 

A transverse view of the warm lead-scintillator calorimeter and its segmentation 
into 88 stacks aligned parallel to the beam pipe can be seen in fig. [2.5| . The stacks 
are multi-layer sandwich structures with active sampling units made of plastic scin- 
tillators. The front face directly behind the BPC is located at z — —144 cm. The 
angular region covered extends from polar angles of 151° up to 176°. 

The electromagnetic energy resolution was derived from test beam measurements 
to be ~ 10%/y/E © 1.7%, the absolute energy scale is known to a precision of 1%. 



'The "©" indicates quadratic addition. 
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Figure 2.4: Central tracking system, sec- 
tion perpendicular to the beam |10 . 



Figure 2.5: Transverse view of the stack 



segmentation of the BEMC [10 . 




Figure 2.6: The orientation of the absorber plates (top) and the cell structure of the LAr 
in (r, z)- view (bottom). 



Chapter 3 

Deep-Inelastic Scattering 



3.1 The Kinematics 

Due to the high beam energies of E e = 27.5 GeV and E p = 820 GeV available at 
HERA, the masses of the electron and proton can safely be neglected for most pur- 
poses. This will be done throughout this study. Their four-momenta,[] indicated in 
fig. pTT| , can therefore be written as k = (E e ,0,0,p z J = (27.5, 0, 0, —27.5) GeV and 
P = (E p ,0,0,p Zp ) = (820,0,0,820) GeV with k 2 = P 2 = 0. Thus, the center of mass 
energy sfs follows from 

s ■= (k + P) 2 = 2k ■ P = 4E e E p = 90200 GeV 2 (3.1) 

to be 300.33 GeV. 

The Neutral Current (NC) interaction displayed in fig. [3.1| is mediated via the 
exchange of a virtual 7 or Z° boson with four-momentum q = k — k' and mass 
q 2 < where k', the four-momentum of the outgoing electron, may be chosen to be 
k' = (E e i, E e i sin 9 e > ,0, E e i cos 9 e >), i.e. 4> e > := 0. Charged Current (CC) reactions involv- 
ing W bosons would yield an outgoing neutrino v e and will not be considered here. In 
the case of elastic scattering ep ep, the Lorentz invariant momentum transfer squared 

Q 2 := -q 2 = -{k - k'f = 2k ■ k' = 2E e E e ,{\ + cos# e /) (3.2) 

would suffice to characterize the process. Protons, however, are not point-like particles, 
but have a complex internal structure revealing itself at distances of < 1 fm, equivalent to 
Q 2 > 0.04 GeV 2 . As consequences the proton, represented by the lower blob in fig. |3.1| , on 
the one hand has to be described in terms of structure functions, and on the other hand 
it usually breaks up since elastic reactions are strongly suppressed compared to inelastic 
ones with increasing Q 2 |T[. This inevitably leads to the need of another quantity to 
define the global outcome of an event, i.e. without differentiating the hadronic final state. 

lr To disentangle four- momenta like p from three- momenta, the latter are set in Roman font with arrows 
on top: p. 
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Figure 3.1: Diagram of the basic ep DIS process via neutral currents. The letters in 
parentheses label the corresponding four-momenta. 



Looking at the photon-proton (7*2?) interaction by itself, we have only one additional 
invariant at our disposal: 

W 2 := (q + P) 2 = 2P-q-Q 2 , (3.3) 

which is limited by < W 2 < s. 

Other popular choices are the two dimensionless variables^ 

* ■■= wi< < 3 - 4 > 

- = U ("> 

with the neat property that < (x,y) < 1. In the context of the infinite momentum 
frame, where the proton is conceived of as a collinear stream of fast moving partons 
and masses are negligible, x can be interpreted as the fraction of the total momentum 
carried by the struck constituent as seen from the virtual boson. This reference frame has 



implicitly been adopted in fig. |3.1| as the hard scattering with a parton, labelled a with 
mass y/~§, is assumed to be incoherent and well separated from ensuing soft processes. In 
the simplest situation of a boson-parton collision where 

s:={q + iP) 2 = 2iq.P-Q 2 = (1 - l\ Q 2 (3.6) 



vanishes, x is identical to £ of fig. EO according to 



l + ^)x. (3.7) 



2 The scaling variable x is also called Bj0rkcn x and y may be denoted as inelasticity. 
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Under more complex circumstances with s > 0, it can only be concluded that x < £ < 1. 

In fixed-target experiments, y is easily interpreted as relative energy loss of the scat- 
tered lepton because in the rest frame of the target denoted by "s: 

y ~ mp E; (3 - 8) 

Of course, only two of the four introduced kinematic quantities are independent of 
each other. The conversion formulae for the pairs (Q 2 , W 2 ) <-> (x,y) are: 



Q 2 Q 2 + W 2 

s 

(3.9) 



X Q 2 + W 2 ' 11 



Q 2 = sxy , W 2 = s{l-x)y. 



3.2 Reconstruction of the Kinematical Quantities 

When confronted with real data, one has to reconstruct the kinematical quantities of an 
event from the measured energy depositions. Basically, we have four measurements at 
our disposal for the determination of two unknowns: the energies and polar angles of the 
scattered electron and the current jet, i.e. a collimated shower of hadrons produced by 
the struck parton: E e r, 9 e >, Ej, 9j. Depending on the choice of input variables, several 
reconstruction methods exist, each with specific advantages and drawbacks. 

Fig. [3.2| presents the lines of constant energy and polar angle for the scattered electron 
and the current jet in the (x, Q 2 )-plane. Lines of constant y = Q 2 /(sx) are displayed for 
y = 1, 10 _1 , 10~ 2 and 10~ 3 . For reconstruction purposes it is best if the isolines of a quan- 
tity are closely staggered and intersect the isolines of the corresponding second variable 
predominantly under large angles. Henceforth, the electron data alone are sufficient for a 
good determination of Q 2 in the complete range shown. The large gaps, however, between 
energies E e > = 25 GeV and E e > = 40 GeV and angles 9 e > larger than 90° demonstrate that 
the electron method is not very reliable for large x or low y respectively. Despite the fact 
that hadronic energies and angles have much larger uncertainties than electromagnetic 
ones, the current jet data can provide better estimates for x or y in that region. 

Yet, because we are interested in the shape of the hadronic final state, we have to 
restrict ourselves mainly to the electron quantities in order not to bias our results. The 



application of a jet algorithm like the ones described in section |4.2.3| to obtain Ej and 9 3 



is not advisable. Fortunately, there is also an inclusive method to derive the kinematics 



from hadronic measurements first proposed by Jacquet and Blondel |14| . Alternatively, 
one could rely on the angular measurements of the electron and the hadronic system 
alone [HJ. These three methods will be presented in the next sections. For a comparison 



ref. [16] may be consulted. 
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Figure 3.2: The two plots show the lines of constant energy (full lines) and polar angle 
(dashed) for the scattered electron (top) and the current jet (bottom) in an (x, Q 2 )-plane. 
Lines of constant y (dotted except for y — 1) are drawn for y = 1, 10 -1 , 1CT 2 and 10 -3 . 
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3.2.1 The Electron Method 

The electron method is the simplest possibility to extract y and Q 2 and employs the 
scattered electron only, x and W 2 can then be calculated according to eqs. (|3.9|) . Using 
the definitions ( |3.2| ) and (|3.5|), it follows directly: 



P-k' Ee>l-cose el E e ,-p z 
y * = l -p-Tk =l -K^^ = l 2E— ' (3 - 10) 



nrp 1 

Ql = 2k-k' = 2E e E e ,{\ + cos0 e /) = 2E e (E e , + p z ,) = . (3.11) 



3.2.2 The Jacquet-Blondel Method 

With h representing the four-momentum of the complete hadronic final state, it is clear 
from fig. |3.1| that q = k — k' = h — P. Furthermore, pt-balance enforces pt h = pt e so that 

P-(h-P) r , 19 . 
Vh= P-k = 2E e ' (3 ' 12) 



2 

Ql = = '- — — (3.13) 

1 - Vh L-Vh 



where % loops over all hadronic objects. 



3.2.3 The Double Angle Method 

Requiring y h = y e and p th = pt , , the dependence on E e > and £\ Ei for all hadrons i can 
be eliminated such that: 

yd = sing e /(l - cosflfe) 

a sin^ e '(l — cos 0^) +sin0 fe (l — cos9 e >) ' 

O 2 - iE 2 sing h (l + cosg e Q 

Vda e sin^(l-cos^)+sin^(l-cos^) ' 1 } 



3.3 The Born Cross Section 

To derive a cross section formula for deep-inelastic ep — > eX scattering, it will be required 
that the total process can be separated into a two step procedure. In the first part, 
involving very small space-time scales of 0(1/Q) where the strong interaction is weak, 
the basic kinematic outcome of a reaction is fixed by the incoherent elastic scattering of 
the boson probe and a proton constituent. The ensuing soft hadronization of the struck 
parton and the proton remnant takes place at a time scale typically of the order of inverse 
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Figure 3.3: Quark-Parton-Model Feynman graph (left) and an example of a virtual cor- 
rection to it (right). 



hadron masses 1/M with M ~ 200 MeV <C Q instead and merely affects the detailed 
structure of the hadronic final state. 

Without resolving such details, it is therefore possible to calculate a cross section by 
applying perturbation theory in lowest or Leading Order (LO, here: O(a 2 a s )) to the 
hard subprocess a of fig. p.lj The blob is then replaced by a single boson-parton vertex 



as shown in the Quark-Parton-Model (QPM) Feynman diagram of fig. |3.3| . If the proton 
itself is left out, essentially an elastic two-body scattering reaction remains for which the 
cross sections have been calculated, s. e.g. One possibility to derive a general cross 
section formula is to assume partons with spin 1/2 and Depending on the spin, one 
takes over the result from e/i respectively e7r scattering: 



2 



dQ 2 Q 4 
da e7T Ana 2 
~dQ^ ~ ~Cf 



;i + (l-2/) 2 ), (3.16) 
[1-y). (3.17) 



Here, a denotes the electromagnetic coupling strength and NC contributions from Z° 
exchange, which are suppressed at least oc Q 2 /(Q 2 + M§), are neglected. 

At this stage, we must again consider the internal arrangement of partons in the proton 
that can be parameterized in the form of parton density functions (pdfs) fi/ p (x) which 
represent the probability to find a constituent % with a momentum fraction in the interval 
[x, x + dx}. Denoting the spin-1/2 pdfs with fi/ p (x) and the spin-0 ones with fj/ p (x), the 
differential cross section can be written as 

= "i 4 " ( 1 + ^2 ~ ^ Qifi/ P (x)dx + (l-y)^2 Qjfj/ P ( x ) dx j , ( 3 - 18 ) 
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where qi and qj are the corresponding electromagnetic charges of the constituent flavours 
i and j in units of the positron charge. Historically, the derivation of the cross section 
implied the evaluation of a leptonic tensor L^ u prescribed by QED and its contraction 
with the most general form of a hadronic tensor W^ u describing the hadronic vertex. This 
approach, employed e.g. in || U, led to the definition of structure functions Fi(x) and 
F 2 (x) which here translate to 

F i( x ) = E^/ P (^ (3-19) 

i 

f 2 (x) = j2^ x M x ) + J2^ x ^/p^- ( 3 - 2 °) 



The double differential cross section Q3.18| ) now reads 

d 2 a ep Ana 2 



dQ 2 dx xQ 4 



((l-y)F 2 (x)+xy 2 F l (x)) . (3.21) 



So far the partons served as a term for point-like constituents within the proton that 
lead to the experimentally observed 1/Q 4 - dependence of the inelastic cross section in 
contrast to the dramatic l/Q 12 -decrease for elastic reactions. Combining this with the 
quark model by identifying partons and spin-1/2 quarks, one immediately concludes that 
fj/ p (x) = OVj and henceforth 

F 2 (x) = 2xF 1 (x) . (3.22) 

This is the experimentally well established Callan-Gross relation which demonstrates the 
fermionic nature of the charged proton constituents. In addition, one can conclude from 
eq. ( |3.21|) that the cross section normalized to the corresponding one for point-like particles 
depends on the scaling variable x only. The observation of deviations can be attributed 
to the neglect of masses, intrinsic transverse momenta and — most importantly — the 
strong interaction. 



3.4 The NLO Cross Section 

Up to now, the proton was treated like a stream of collinear non-interacting quarks. 
However, in Next-to-Leading Order (NLO), i.e. 0(a 2 al), quarks can emit and absorb 
gluons which again may split up into two gluons or qq pairs and so forth. The rather 
simple picture involving structure functions Fi(x), F 2 (x) to describe the composition of 
the proton in terms of quark densities alone has to be modified accordingly to include 
a gluon density. Two of the most important consequences are scaling violations, i.e. 
a dependence of Fi(x),F 2 (x) on Q 2 , and non-vanishing contributions to the longitudinal 
structure function 

F L := F 2 - 2xF 1 . (3.23) 
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Figure 3.4: QCD-Compton Feynman graphs. 



For low x, i.e. x < 10 3 , approximations reveal the effects to be proportional to a s (Q 2 
times the gluon density g(x, Q 2 ): 



dF{x, Q 2 
dlogQ 2 



oc a s (Q 2 ) ■ g(x,Q 2 



(3.24) 



F L (x,Q 2 )^a s (Q 2 )-g(x,Q 2 ). (3.25) 

Both are intensively studied in DIS experiments and can be exploited in inclusive mea- 
surements to gain information on the strong coupling constant as well as the gluon density. 
For HI publications on this topic consult refs. JT7], P~8| . 

In this analysis all hadronic final states are considered but will be differentiated with 
respect to suitably chosen characteristic variables (s. section [4.2[ ), i.e. the measurement is 
semi-inclusive. 

To 0(a).) there are two kinds of processes resulting in two instead of one final state 
parton, where "parton" from now on is used synonymously for quarks, anti-quarks and 
gluons. Since they may lead to "real" effects like an additional jet they are called real 
corrections. In the first case of the QCD-Compton graphs (QCDC), an additional gluon 
is emitted by the struck quark. The two possible Feynman diagrams are presented in 

involves the production of a qq 



fig. |3^j. The second kind with diagrams shown in fig. |3T 
pair from a reaction between the boson and a gluon emitted from the proton. They are 
expressively labelled Boson-Gluon-Fusion graphs (BGF). 

To complete the set of contributions to ep DIS to 0(a]), one also has to consider loop 
diagrams. One example is given on the right-hand side of fig. [3.3| . The diagram itself 
results in an 0(a 2 ) add-on to the matrix element squared. However, since the final states 
of both graphs in fig. |3.3| are indistinguishable, an interference term of 0(a]) has to be 
taken into account. For the purpose of differentiating the hadronic final state, this virtual 
correction is, of course, irrelevant and merely changes the total cross section. 
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Figure 3.5: Boson-Gluon-Fusion Feynman graphs. 



Focusing on the two-parton QCDC and BGF processes, three new degrees of freedom 
arise in the matrix elements corresponding to energy, azimuthal and polar angle of one 
parton. The other parton is then determined by energy-momentum conservation. Usually, 
the now five-fold differential cross section is written employing the variables 



Zr := £_^ =1 -£_£5, (3.27) 
F ■ q P ■ q 

and (f>, where <fi denotes the angle between the planes fixed by the two outgoing partons on 
the one hand, and the scattered electron on the other hand in a suitable reference frame. 
This may be e.g. the Breit frame defined in the next chapter or any other reference system 
connected to it by a Lorentz boost along its z-direction. Concerning the characterization 
of the hadronic final state, no reference will be made to the electron. Therefore, the 
integration over can be performed beforehand such that for our purposes the NLO cross 
section may be labelled as 

d 4 a ep 

(3.28) 



The allowed ranges for x p and z v are x < x p < 1 and < z p < 1. 

Yet, it does contain singularities! For the QCDC and BGF processes they can be 
extracted from ref. [T9], which presents the complete 0{a 1 s ) corrections to electroweak 



NC and CC ep DIS cross sections, to be: 
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d " QCDC « (3 ' 29) 

da BGF oc ^ + (1 " { (1 ~ ^ p)2] . (3.30) 

Partially, they cancel against corresponding divergent terms of the virtual corrections. 
The remaining initial state mass singularities can be absorbed in a proper redefinition 
of the parton densities. Nevertheless, one has to be careful about the precise definition 
of quantities that shall be calculated by pQCD. In order to allow the divergences to 
compensate each other, the investigated variables must agree in the limits where the real 
and virtual corrections go to infinity. The number of final state particles for example 
would always be two in the first case and one in the latter. 

To be more definite, we make use of eqs. ( |4.22j )— fl4.26|) and can refer the limits x p — > 1 



and z p — > 0, 1 to configurations where either p\ or p2 are soft (infrared divergence) or any 
pair of {pi,p2,r} is collinear (collinear divergence). Hence, any quantity F n (pi, . . . , p n ) 
to be meaningful in pQCD has to fulfil the conditions 

F n (pi,..., Api, ...,j>j = (1 - A)pi, . . . ,p n ) = Fn-lCpi, • • • ,Pi, • • • ,Pj-l,Pj+l, • • • ,p n ) 

and (3.31) 
limF„(pi,...,Api,...,p n ) = F n _i(pi,...,p i _ 1 ,p i+ i,...,p n ) 

A — ^0 

for < A < 1. That is, collinear splittings and soft particles do not affect F; it is collinear 
and infrared safe. One example for such a variable is thrust, or rather 1— thrust, defined 
by eq. ( f4.9| ). It is invariant with respect to collinear splittings and varies smoothly for one 
momentum approaching zero. 



Chapter 4 
Event Shapes 



4.1 The Breit Frame 



4.1.1 Introduction 

Initially, all energy deposits measured with the HI detector are given in the laboratory 
system. Since we are, however, merely interested in specific properties of the hadronic 
final state caused by the hard interaction, two problems arise: 



1. The transverse momentum of the scattered electron is balanced by the hadronic 
system. Hence, from longitudinal and transverse momenta of the hadrons one can 
reconstruct the global event kinematics, but they are not characteristic of the un- 
derlying hadronic process. 



2. Somehow we have to differentiate between the products of the hard reaction and 
the proton remnant. A maximal separation, which is not given in the laboratory 
system, is desirable. 



For illustration fig. [O] shows a comparatively simple NC event with a clearly identifiable 
electron and a lot of hadronic activity on the opposite side with respect to <fi. In addition, 
the proton rest manifests itself in the form of some clusters in the forward direction. 

Here already, it is not too obvious how to define characteristic properties of the ha- 
dronic energy deposits. In the case of much more complicated average events, it becomes 
forbidding. The solution to the two problems is to apply a Lorentz transformation into the 
Breit frame of reference [2~0f|. For an event of QPM type as shown in fig. |4.2| , it is defined 



as the reference system where the incoming parton with momentum xP* in +2;*-direction[] 
is back-scattered by a purely space-like boson of momentum 



1 Note that this is in contrast to e.g. (|o) and [^lj where the +z*-axis has been chosen for q*. 

2 To distinguish non-invariant quantities in the Breit system from those in the laboratory, they will be 



marked by a *. 
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Figure 4.1: Example of an ep collision measured with the HI detector which resembles a 
QPM-like configuration. The isolated electromagnetic cluster in the LAr with one track 
linked to it represents the scattered electron, whereas the broader energy deposits in the 
upper half of the LAr balance the electron p t and correspond to a jet of hadrons produced 
by the struck quark. Some energy measured in the very forward direction indicates the 
proton remnant. The basic kinematic quantities deduced from the electron are given. 

From g* 2 = q 2 = —Q 2 and xP* + q* = —xP*, it follows that 

g* = (0,0,0,-Q) (4.1) 

and 

where Q is defined to be 

Q := \Tq~ 2 = \[Q 2 . (4-3) 

Thereby, the event and also the "soft" and "hard" physics is separated by the (x*,y*)- 
plane into a +z* Remnant Hemisphere (RH) and a —z* Current Hemisphere (CH). 
Simultaneously, the transverse momentum of the scattered quark has been eliminated such 
that the CH with its available energy of Q/2 is very similar to one half of an e + e~ — * qq 
event with purely "time-like" energy y/s/2, relating the "relevant energies" Q and y/s. 
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Figure 4.2: Diagram of a QPM-type ep collision as seen in the Breit frame. 



Technically, the necessary Lorentz transformation is decomposed into a pure boost 
demanding that 

2xP* + q* = ( ? ] (4.4) 



and a rotation afterwards to readjust P* into the +z*-direction. In addition, the incoming 
and outgoing electron is usually required to lie in the (x*, z*)-plane. Their boosted and 
rotated four-momenta read 



and 



*. = l |(2-»)«^,0,-f 



(4.5) 



(4.6) 



Thus, another feature of the Breit frame is that the energy loss of the lepton, E* — E e *', 
vanishes. 



4.1.2 Properties 

For a better clarification of what happens to the clusters in fig. [4.1| , we introduce another 
variable 9 q which corresponds to the polar angle of the scattered quark as seen by the 
electron: 

cos 9 = 910^ZM3 (47) 

A comparison of 9 q , drawn also in fig. |4.1| , with the polar angle of the most energetic ha- 
dronic energy deposition in the LAr obviously demonstrates that they are approximately 
equal. When looked at it from the Breit frame, the +z-axis remains at its position, but 
the new — z*-direction is given by 9 q \ As a result, the angular region of 0° < 9 < 9 q is 
stretched until 0° < 9* < 180° and, together with the remaining 180° — 9 q , the complete 
polar angular range for opposite 0, including the electron and most of the proton remnant, 
is squeezed. Intermediate angles other from (fi q and e with (fi q — e = 180° constitute 
the transitional domain between the two cases. 

For all variables defined in the next section, the energies and polar angles of the trans- 
formed four-vectors are the important input quantities. Azimuthal angles with respect 



to that of the scattered electron, i.e. 0* = 0°, are rather insignificant. In fig. 4.3 we 
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therefore restrict ourselves to four plots showing the 6*(9,(j))-, 8*(9,E)-, E*(E,<p)- and 
E*(E, #)-functions for the sample boost of fig. 4.1 . For simplicity, an exact balance in p t 
of the electron and the hadronic final state is assumed so that E = 85.5 GeV is taken as 
default. 

From these plots it can be concluded that: 



1. 9*(9,4>) peaks sharply at 



'qi YqJ 



2. Only clusters opposite of the electron hemisphere in the laboratory system appear 
in the CH, i.e. 9* > 90°. 

3. 9* is independent of the energy of the boosted vector. 

4. E* oc E with the slope strongest at e and 9 = 180°. 

The first property may lead to a serious deterioration of the resolution in polar angle in 
the Breit frame and is the motivation for the cut-off no. |5| introduced in section |5.5| . 



4.2 Definition of the Event Shapes 



According to fig. |4.2j , the hadronic final state in the simplest case consists of merely 
one parton with longitudinal momentum only and no mass. The hard interaction is of 
a purely electromagnetic nature. QCD induces deviations from this constellation. By 
investigating variables that are sensitive to these deviations, it is possible to learn more 
about perturbative and non-perturbative aspects of QCD. 

All quantities, introduced in the next sections and generically labelled as F, are event 
shapes and have the following properties in common: 

1. They are dimensionless. 

2. For ep DIS they are defined in the Breit frame of reference. 

3. In the limit of a QPM-like event F — 0, otherwise F > 0. 

4. They are infrared and collinear safe, which is most important for a valid comparison 
with pQCD. 

5. Soft fragmentation and hadronization processes cause discrepancies between fixed- 
order pQCD calculations and measured data that generally can be parameterized 
to be oc 1/Q P with Q being the relevant energy scale and p = 1 or p = 2 in our 
circumstances. 

The event shapes discussed in the next two subsections assume that the hard sub- 
process takes place in the CH alone. However, for £ > 2x it is possible for the CH to 
be completely empty, although experimentally this is improbable due to hadronization, 
backscattering, noise, etc. In order to be insensitive to such effects and to remain infrared 




Figure 4.3: The 9*(0,(f>)- (top left), 9*(6,E)- (top right), E*(E,cf>)- (bottom left) and 
E*(E, #)-functions (bottom right) for the sample boost of fig. |4.1| . The third input variable 
is set to E — 85.5 GeV, 4> = (fi q = 111.3°, 6 = 9 q = 37.3° and again (ft = (f> q respectively. 
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safe P2| , one has to specify what is meant by an "empty" or, vice versa, a "full" CH. We 
adopt the prescription that the total energy available in the CH E* := ^2 iecn E* has to 
exceed 20% of the value it should have according to QPM: 

E* > 0.2 • Q/2 = Q/10 . (4.8) 

Otherwise the event is ignored. The precise value of the cut-off is motivated by a study 



of the measured energy flow in both hemispheres performed for ref. [^3| . 

To keep event shapes dimensionless, it was originally suggested in [[24] to normalize en- 
ergies, momenta and masses to Q/2. "Empty" events would then lead to F — 0. Following 
the proposal in |25[ to use for experimental reasons the actually present total energy E* 
or momentum P* := E; £ ch |pV| instead, one would get ill-defined expressions. Therefore, 
F is set to zero in [[H] analogously for this normalization in contrast to our definition, 
thereby affecting the total cross section a and the left-most bin of the differential cross 
section da/dF. Except for the difference in a, mean values (F), however, are not altered. 
In this work, we will mainly restrict ourselves to the study of the event shapes normalized 
to E* or P*. 

Note that in contrast to the first published experimental results on event shapes in ep 



DIS [p25[ f23f we adopt the modified naming scheme from ref. ||21|| . Except for tc defined 
below, the subscript indicates the quantity (E, P or Q) used for the normalization. The 
event shapes 1 — T c , 1 — T z , B c and pc investigated in the above-mentioned publications 
will be labelled t c , t p , B p and pn. 



4.2.1 Event Shapes employing z* as Event Axis 

Choosing the direction of the exchanged boson q * = (0, 0, —Q) as event axis n, one can, 



with fig. [4.2| in mind, easily deduce quantities F being zero for QPM-like configurations 
and F > otherwise. The simplest event shapes thrust (or rather 1— thrust r := 1 — T) 
and jet broadening B can be written as [|(| |24], [25] 

E |pi*-s| £ \pu\ E \pu\ 

Tp:=1 "^wr =1 "~^' tq ^ x --qit> (49) 

iSCH 

E W x n| E \Pti\ E \Pu\ 

p := 2 e m = ' := ~q~ (4 - 10) 

igCH 



where pi and pt denote the longitudinal respectively the transversal momentum compo- 
nents of p. The factor of 1/2 for B is conventional. 
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4.2.2 Event Shapes without Reference to z* as Event Axis 

A possibility for differentiation without reference to z* as event axis is given by the jet 
mass p p7| : 



Pe 



;sch 



M 2 



4 ( E Ei 



AE 



c2 • 



M 2 



with M being the total mass in the CH. 

Another way is the evaluation with respect to a new axis to be defined. Originally, 
event shapes were introduced in the context of e + e~ annihilations, where, in distinction 
to ep DIS, a preferred direction like z* is not given beforehand. The first definition deals 
with the thrust axis nV which is characterized as the normalized vector that maximizes 
the sum of all projections of momenta (absolute values) onto it. With regard to that, 
t = 1— thrust can now be described as 



E |Pi*-n| E |Pi*-n T | 

t c := 1 - max ' 6C " - , = 1 - . ( 4 .12) 

a,tf»=i E Pi* P* 



iGCH 



In contrast to e + e physics, the maximizing procedure is applied in the CH alone. For 
merely one momentum vector, Tc always equals zero. 

A second kind of event axis brings the momentum tensor 



igCH 1^ 

ieCH 



(4.13) 



into play. This real symmetric matrix is positive semi-definite with trace Tr(@) = Ai + 
A2 + A3 = 1. For < A3 < A2 < Ai < 1 it describes an ellipsoid with pairwise orthogonal 
axes named minor, semi-major and major with increasing eigenvalue. The major axis is 
similar but not identical to nV. If A3 = 0, then the ellipsoid degenerates into an elliptical 
cylinder with all momentum vectors in one plane. Is A2 = as well, then all momenta 
are collinear and the corresponding normal area consists of one or two parallel planes. 
Utilizing the eigenvalues, we can define the C parameter [p£j[| : 



C P := 3(AxA 2 + A 2 A 3 + A 3 Ax) , C Q := . (4.14) 



The conventional factor of three ensures a maximal value of one for Cp. 
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4.2.3 Event Shapes employing Jet Algorithms for the Separa- 
tion of the Remnant 

Another approach of characterizing an event with regard to deviations from the QPM type 



does not make use of the CH. As depicted in section |3.1] , the elementary reaction yields 
one hard parton moving along the — z*-direction. The ensuing production of soft partons 
during fragmentation and the final hadronization distort this picture to a limited extent 
only since no large transverse momenta with respect to the original one are involved. 
Basically, the one parton gets transformed into a tight stream of hadrons, flying along the 
original direction, which one refers to as a jet. 

The inclusion of more complex processes into the pQCD calculation facilitates the 



production of two or more hard partons (s. section 3^), and henceforth the events can 
acquire more than one jet in addition to the remnant jet. QPM-like events may also be 
called to be of a (1 + l)-jet type in contrast to (n + l)-jet events with n > 2. In order 
to decide to which category a given constellation belongs, precise instructions on how to 
combine jets from an assortment of four-momenta are needed. The most basic of these jet 



algorithms [£9[ makes use of angular cones around seeds given by the input four-vectors. 
The two schemes we employ are of another type called cluster algorithms. Both, the 
Durham- or k t - [|30 , |3i| and the JADE-algorithm [p2j [331 ar e applied in a modified form 



adapted to ep DIS in the Breit frame. 

The central procedure is almost the same for both. Two distance measures are defined, 
one for distances between two four- vectors, yy, and another one for the separation of each 
from the remnant, yi r . When all combinations are evaluated, the minimal value determines 
which are the closest two in y. If an was smallest, then these two are recombined to 
one new four- vector. In case of yi r to be minimal, % is ascribed to the remnant. The 
whole routine is repeated until either all y's are larger than a lower bound y cu t, or until 
a certain number of jets is reached. The first prescription is used to divide a sample of 
reactions into sets of (1 + 1)-, (2 + 1)- and so forth events. The second approach is taken 
to employ y as an event shape variable. Here, y^ t and yjj always denote the y- value where 
the transition (2 + 1) — > (1 + 1) occurs. The respective distance measures are 

Vii ■= 3 Q2 -> (4-15) 



2^^(1 -008 6?; 



Q 2 



(4.16) 



for the factorizable JADE-algorithm and 



2min(^,^)(l-cos^ v 

Vij ■= Q- 2 ' ( 4 - 17 ) 



for the Durham-algorithm. 
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4.3 Event Shapes to 0{ct s ) 

As explained above, the event shapes are designed to distinguish between a QPM-like 
topology of the hadronic final state with the relevant cross section given by eq. ( 3.21 ) and 
deviations from it that are described in lowest order pQCD by eq. ( 3.28|) . In the first case, 
F always equals zero, whereas in the latter, F depends on the two additional degrees of 
freedom x p and z p . 

To derive explicit formulae, we first determine p\ and P2 from 



Pi+P 2 = £ p +1 



( SL \ 

o 





(4.19) 



according to eqs. ( |4.1| ) and ( |4.2| ). Remembering that azimuthal angles are irrelevant, we 
define the y-components to be zero such that 



Using p\ 
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(4.20) 
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z — z 3 



Zq — Z 3 



one can calculate Zq, z , z t , z 3 and z 3 to be 



Z() 

Zq 
Zt 
Z3 
Z3 



2,XpZp 



Xt 



z v + l 



,1, J 



'^CCpZp ~\~ Xp ~\~ Zp 



X 



2j^(l- Xp )(l-z p ) 

Xp 

Xn Zn 



Zp Xp 

X r . 



In the limit of x r , — > 1 and z„ 



(4.21) 



(4.22) 
(4.23) 
(4.24) 
(4.25) 
(4.26) 





1 they evaluate to zq = 1, z 3 = — 1 and zq = z 3 = z t 
corresponding to a QPM-like event with only one final parton. Taking into account the 
mismatch in the sign of the z*-direction and noting that x p is named £ in [^TJ, the formulae 
coincide with those of [20| and [f2l|1 . 
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Figure 4.4: (x p , z p ) phase space region for 0(a s ) corrections to DIS. See text for an 
explanation of the labels. 

Finally, to achieve results for the event shapes, the available phase space in (x p , z p ) 
has to be subdivided. Fig. [4.4| shows on the left-hand side the appropriate subregions for 



the definitions of sections f4.2.1| and [4.2.2| involving the CH only. The four triangles A to 
D correspond to: 

A: Both partons are in the CH: z 3 ,z 3 < 0. 
B: Only parton p* is in the CH: z 3 < 0, z 3 > 0. 
C: Only parton p\ ls m the CH: z 3 > 0, z 3 < 0. 
D: The CH is empty: z 3 , z 3 > 0. 

Except for subregion D which is excluded from our event shape definitions by the 
explicit requirement of a minimal energy of Q/10 in the CH, all results are listed in 



table |4.1| . The dotted lines in fig. fl4| point out two parts of regions B and C rejected in 
addition to D. Note that p#, pg, Cp and Cq vanish throughout B and C. 

The additional separation into Ai and A2 indicated by the dashed line is necessary for 
t c . Depending on the angle enclosed by pi* and pV, the thrust axis n T is represented by 



n^ 



pi*+iV 

Pl*+j32* 



Pi ~P2 
|Pl*-P2* 



for Z 12 < 90° (Ai) and 
for Z12 > 90° (A 2 ) . 



(4.27) 



In the first case, tq is equal to rp, in the latter, tc also measures momentum components 
perpendicular to the boson axis n as can be seen from the appearance of Zt in the formula 
for A 2 . 



For the event shapes employing jet algorithms as explained in section |4.2.3| , the com 



plete phase space is accessed. As displayed on the right-hand side of fig. [14], it is split up 
into three main regions: 
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A: The two partons are merged together. 
B: Parton p\ is clustered to the remnant. 
C: Parton p\ is clustered to the remnant. 

The full line designates the border between A, B and C for y^, whereas the dashed line is 
valid for yjj. Here, the subdivision into Aj and A 2 reflects the min(.Ej* 2 , Ej* 2 ) condition 
of eq. (pTTTD . 



It is remarkable that to 0{a s ) some of the definitions from section fO| lead to the 
same formulae, e.g. for tq, pg and yjj. However, considering the complete phase space 
in (x p ,Zp), discrepancies appear. When higher orders are included, all event shapes will 



differ from each other. Table fO| gives an overview of the allowed ranges for the defined 
event shapes. 
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F 
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T Q 
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z t /2/z 
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Zt 


z t /2 
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PQ 


(1 Xp) / Xp 
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2(1 -Xp) 1 - ^4xjz 2 + (1 - 2z p f 








c P 


3(2x p - l) 2 z 2 /(z z ) 








Cq 


3(2x p - l) 2 z 2 /(x 2 p z z ) 








VfJ 


(1 Xp ) 1 Xp 


Zp 


(1 - Zp) 


Vkt 


z /z ■ (1 - Xp)/X p Zq/Zq ■ (1 - Xp)/X p 


ZpZo 


(1 - z p )z 



Table 4.1: Formulae for the event shapes F in 0(a s ). For the definition of the phase 
space subdivisions A 1; A 2 , B and C see fig. fO| and the explanations in the text. 



F 


Upper bounds 
to 0(a s ) absolute 


F 


Upper bounds 
to 0(a s ) absolute 


T~P 


1 


1 


T Q 


1 


1 


Bp 


1/2 


1/2 


Bq 


1 


1* 


PE 


1/4 


1/4 


PQ 


1 


1* 




1 - V2/2 


1/2 








c P 


3/4 


1 


Cq 


3( v / 5-l) 3 (3-v / 5)/4 


2.89* 


VfJ 


1/2 


0.84* 


Vk t 


< l/(4x) 


276* 



Table 4.2: Upper bounds for the defined event shapes. The starred numbers are not 
absolute values, but the largest ones encountered. 



Chapter 5 
Data Selection 



The very basis of every collider experiment is formed by the equation 

dN „ . . 

- = £.. (5.1) 

It relates the observed event rate dN/dt with the corresponding cross section a. The 
machine dependent proportional factor £ is called luminosity and has to be measured 
e.g. via comparison to a theoretically well-known reaction like the Bethe-Heitler process 
ep — > ep7 employed in HI (s. section |2.2.1|) . 

In order to gather as many events of a certain kind as possible, one would like to have 
a large luminosity. For Gaussian beam profiles with horizontal and vertical widths o~* , 
°* yi and a* 2 , it is given by 

£ = 4t—s — V (5-2) 

where 1%, I2 are the beam currents and /& is the bunch frequency. By increasing the 
currents, the HERA crew was able to improve the performance considerably over the 



years. Fig. |5J] gives an overview of the integrated luminosity £; n t that was produced and 



accumulated during the running periods from 1992 up to 1997.[| 



5.1 Background Sources 

Nevertheless, it should be kept in mind that background processes are enhanced right 
along. In fact, interactions with atoms of the rest gas in the beam pipe are even dominat- 
ing! Table gives an impression of the rates to expect. Therefore, a careful consideration 
of possible background sources is necessary. Basically, they can be subdivided into "true" 
background not related to ep collisions and misinterpreted competing ep reactions. The 
first are common to all analyses, whereas the latter depend on the topic under study. 

1 Cross sections a are usually measured in barns b, where 1 b = 10~ 24 cm~ 2 . The luminosity integrated 
over the time £j nt may therefore be given in e.g. pb _1 . 



37 



38 



CHAPTER 5. DATA SELECTION 




Cross sections and rates (at design luminosity) 


Beam gas interactions 




50000 Hz 


Cosmic myons 




700 Hz 


Photoproduction 


1600 nb 


25 Hz 


NC DIS, Q 2 < 100 GeV 2 


150 nb 


2.2 Hz 


NC DIS, Q 2 > 100 GeV 2 


1.5 nb 


0.022 Hz 


CC DIS, P t > 25 GeV 


0.05 nb 


0.001 Hz 



Table 5.1: Cross sections and rates at a design luminosity of 1.5 x 10 31 cm 2 s 1 [|J. 
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5.1.1 Non-ep Background 

Beam Gas Interactions 

Instead of colliding with particles from the other beam, it is also possible (and even 
probable) to hit residual gas atoms. Especially the high energetic protons are able to 
produce large numbers of particles that may be scattered into the detector. Most of these 
events can be rejected because there is no vertex or the tracks point to vertices outside 
the interaction region. 



Cosmic Myons 



The surface of the earth and henceforth the detector are constantly hit by myons of cosmic 
origin. Most of them cross the experiment out of time with respect to ep collisions and 
out of place, i.e. nowhere near the interaction region. Owing to their high rate, however, 
sometimes they pass in time right through the CTD and fake DIS events, despite the 
fact that they usually deposit little energy in the LAr. Still, they can be identified 
topologically by looking for back to back tracks in the IRON and CTD. In addition, the 
produced clusters are low energetic and imbalanced in transverse momentum. 



Halo Myons 

A second source of myons are stray protons interacting with material around the beam line. 
They always surround the proton beam (halo) and occasionally give rise to electromagnetic 
showers in the calorimeters, but since they are moving along the proton direction, they 
usually do not cause high p t tracks in the CJC. 



5.1.2 ep Background 

Concerning our aim to investigate NC DIS, the following processes have to be considered 
as background. 



Photoproduction 

The term "photoproduction" derives from the picture of a quasi-real photon interacting 
with the proton, i.e. Q 2 ~ 0. Hence, the electron is only slightly deflected from its 
original direction and escapes the central detectors. Instead, it can be found with a 
certain efficiency in the electron tagger of the luminosity system. Nevertheless, the ^p 
reaction may result in a hadronic final state involving large transverse momenta. The 
misidentification of a particle as the scattered electron is then, due to the very high rate 



s. table 5.1), a main source of background. 
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Diffractive Events 

Another name for these reactions is "Large Rapidity Gap events" (LRG's), where the 
rapidity is defined as 

1 ^ E + P* 

yR:= 2 ln E^- z - (5 ' 3) 



In the case of massless four- moment a, it is identical to the pseudo-rapidity 

7} := - In ^ tan • ( 5 - 4 ) 

In the region between the proton remnant and the current jet, which are normally con- 
nected via colour forces, an unusual gap without hadronic activity is exhibited. These 
events constitute around 5-10% of the NC DIS sample, yet they are not described by the 
pQCD calculations invoked and have to be excluded. 



CC Reactions 

At last, it may happen that clusters of the hadronic final state in one of the seldom CC 
events with a neutrino as scattered lepton are misidentified as electron. Since the neutrino, 
however, escapes detection completely, the total transverse momentum is imbalanced. 



5.2 Trigger Scheme 

To keep notwithstanding the high rates the dead time of the experiment low and specif- 
ically select true ep collisions, a trigger system has been set up in four stages ||. The 
total read-out time of the detector of ~ 1 ms is four orders of magnitude larger than that 
between two bunch-crossings of ~ 0.1/xs. Therefore, only a small part of the measured 
data, especially those from detector parts with short response times, is available for a fast 
decision. 

The first level trigger LI collects information from nine trigger systems attached to 
one subdetector each. These trigger elements are combined to form various subtriggers 
which provide a KEEP or REJECT signal within 2.5 /is. A pipelining system stores the 
full data at the front end during the delay caused by LI and ensures a dead time free 
running at this stage. 

For future requirements intermediate trigger levels L2(] and L3 are foreseen to operate 
during primary dead time of the read-out and are based on the same input as LI. However, 
they are able to evaluate a larger number of signals and their complex correlations. The 
decision times of these systems are designed to be around 20 /xs and 800 /is respectively. 

The last stage consists of the L4 software trigger, which has the raw data of the full 
event at its disposal. Depending on the time consumption, either fast algorithms specif- 
ically adapted to the requirements of L4 or parts of the standard offline reconstruction 

2 L2 was commissioned in 1996. 
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program H1REC []33J are applied to reach a quick decision. All events accepted by L4 are 
finally stored on tape. In addition, a small fraction of L4 ReJecTed events (L4RJT's) of 
~ 1% are kept for monitoring purposes. 

The complete scheme including a fifth level L5 discussed in the next section is illus- 
trated in fig. [57|. 

5.3 Event Classification 

Even after all four trigger stages, the reconstructed events kept for physics analysis contain 
to a large part unwanted or background reactions for most investigations. To save tape, 
network and computing resources, the program FPACK |35| used for platform independent 
data access has the possibility to skip events according to a classification word stored 
directly at the beginning of an event. Thereby, only the properly marked data tracks are 
completely read. Depending on the aim of a study, different classes, each matching to a set 
bit in the classification word, can be selected. Events falling in no physics class at all are 
L5 ReJecTed (L5RJT). They are kept as raw data, but do not appear fully reconstructed 
on Production Output Tapes (POT's), except for a small fraction for monitoring purposes 
as in the case of L4RJT's. 

The chief requirement distinguishing NC DIS from background processes is the pres- 
ence of a candidate for the scattered electron in the BEMC or LAr calorimeter. Therefore, 
loose criteria such as the existence of at least one compact mainly electromagnetic cluster 
with a minimal energy of several GeV are applied to define two DIS classes, a low Q 2 
class (no. 11) in case of the BEMC and a high Q 2 one (no. 9) for the LAr. 

Using additionally data from the tracking and myon systems, beam-induced and cos- 
mic background is rejected more effectively. 

5.4 Preselection 

Due to the soft cuts applied in the event classification, the resulting classes are usually 
still too contaminated with background. In addition, some analyses might not need all 
parts of the intended content. Yet, physics studies are refined and iterated several times 
before being finished so it is recommendable to reduce the data sample further. This is 
done by a preselection which copies the selected events to a local disk for fast access. 

Since the operation of such a complicated machinery like HERA and HI is no easy task, 
running conditions are subject to variations. Any time a significant change in status like 
the (temporary) failure of an important detector component occurs, a new run is initiated. 
Only runs qualified as "medium" or "good" are allowed by the preselection. Additionally, 
the most important detector components for this analysis, i.e. the CJC1/2 and LAr for 
class 9 and also the BEMC and BPC for class 11, are required to be operational. 

Relying on the fact that at least the scattered electron should cause a well measured 
track, the existence of a reconstructed vertex is mandatory. 
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Figure 5.2: A principal diagram of the data stream in HI. Three levels of triggers operat- 
ing with fast signals from nine detector branches and one software trigger evaluating the 
complete raw data are responsible for the rejection of background. A fifth level L5 works 
on fully reconstructed events and categorizes them into physics classes that are stored as 
POT's. To avoid overflowing them, the classes may be downscaled. Events accepted by 
L4 but falling in no physics class at all are kept as raw data only. Small samples of L4 
and L5 rejected data are retained for monitoring purposes. 
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As mentioned above, it is evident for NC DIS that finally a clear candidate for the 
scattered electron should be found. In distinction to the event classification, stricter cuts 
have to be fulfilled now. In HI a program package for PHysics ANalysis (H1PHAN PS| ) 
is available containing several algorithms for that purpose. Two of them, QFSELH/M 
and QESCAT, are applied. Both exploit the properties that the scattered electron should 
be isolated from the jet of the struck parton due to pt-balance and that it deposits its 
energy in form of a compact electromagnetic shower mainly in the ECAL part of the LAr 
or in the BEMC. 

In the case of QFSELH/M, some additional requirements are imposed on the found 
candidate. A cone of radius 0.5 in azimuthal angle <fi and pseudo-rapidity 77 is drawn 
around the candidate. At most 5% of its energy is allowed to be measured within this 
cone. Additionally, a matching track has to be found within 5° for high Q 2 . For electrons 
in the BEMC, a hit in the BPC with a maximal distance of 5 cm from the cluster center 
projected onto the BPC (x, y)-plane has to be measured and the cluster radius may not 
exceed 5 cm. 

For the purpose of the preselection, events are accepted if a candidate is determined ei- 
ther way. In the final analysis QFSELH/M is utilized for the low Q 2 sample and QESCAT 
for the high Q 2 events. A comparison revealed marginal differences between the two elec- 
tron finders. 

At last, a minimal energy of the candidate of 9 GeV is called for, and it is ensured 
that safe regions of the calorimeters are hit by cuts on the polar angle of 10° < 9 < 150° 
for the LAr and 157° < 9 < 173° for the BEMC respectively. Because the minimal Q 2 
demanded later is 49 GeV 2 , a slightly smaller value of 45 GeV 2 , calculated according to 
eq. ( |3.11| ), has to be surpassed by the low Q 2 events. 



5.5 Final Cut Scenario 

Before the last selection is applied, an additional routine rejects residual events due to 
cosmic myons. Furthermore, there may be no photons detected in the PD. 

In order to achieve a clean DIS sample, the final cuts have to take into account 
detector acceptances, efficiencies and resolutions as well as the diverse background sources 
described in section The latter have already been suppressed to some extent in the 



selection procedure above. Additionally, possible constraints of the theoretical model to 
compare with have to be considered. 

The final cuts are subdivided into two basic categories: 

• Phase space cuts: 

Due to unavoidable limitations either of the experimental apparatus or the theoret- 
ical model, it is not possible to choose at will a phase space to investigate. Hence, 
common kinematic requirements have to be imposed on the data as well as on the 
theory. 
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Data quality cuts: 

These are necessary to ensure that background is suppressed and the measurement 
is of good quality. Then, the data can be corrected for detector effects and may be 
compared to theoretical predictions. 



First, the phase space cuts will be defined. For a quick overview see table |5\2 . 
Phase space cuts: 

• low Q 2 : 49 < Q 2 / GeV 2 < 100, high Q 2 : 196 < Q 2 / GeV 2 < 10000: 

To examine the Q-dependence of the event shapes, all events are grouped into eight 
bins in Q: 7-8 GeV, 8-10 GeV, 14-16 GeV, 16-20 GeV, 20-30 GeV, 30-50 GeV, 
50-70 GeV and 70-100 GeV. Note that there is a gap from 10 GeV to 14 GeV 
which corresponds to the excluded transition region between the BEMC and LAr 
calorimeters. The lower bound of 7 GeV is motivated by the fact that an energy of 
Q/2 = 3.5 GeV should be available in the CH of the Breit frame. For Q > 100 GeV 
there is not enough statistics at hand. 

In case of data, Q 2 is identified with Q\ according to eq. (|3.11|) . 



0.05 <y < 0.8: 

As explained in section |3.2| , the electron method is not well suited for the recon- 
struction of y-values as low as 0.05. Hence, these events are rejected. At high y, 
radiative corrections due to an additional real photon emission of the scattered elec- 
tron (s. section |6.1.1| ) are enormous and have to be avoided. 

For data, y is identified with y e according to eq. ( |3.1U| ) down to 0.15. Below that 



value, yh from eq. (|3.12|) is taken instead.^ Note that nevertheless y e > 0.05 is asked 
for. 

low Q 2 : E e , > 14 GeV, high Q 2 : E e , > 11 GeV: 

Lower limits on the electron energy are imposed for several reasons: First, the trigger 



efficiency is above 99% [37, [13| . Moreover, there are numerous clusters of hadronic 
origin, e.g. pions, which are low energetic and may fake electrons, especially when 
the real one did not even hit the main detectors e.g. in photoproduction. At last, 
true leptons which lost a large part of their original energy due to radiative effects 
are excluded. Thereby, a good measurement of E e > should be reached, which is most 
important for the boost into the Breit frame. 

low Q 2 : 157° < e , < 173°, high Q 2 : 30° < e , < 150°: 

These cuts reflect the coverage in polar angle of the BEMC and LAr calorime- 
ters, although the upper limit of 173° for the low Q 2 sample is redundant owing to 
Q 2 > 49 GeV 2 . For the high Q 2 events, the forward region 9 e / < 30° with its high 
hadronic activity is left out to avoid misidentifications of the scattered electron. 



^Concerning the boost into the Breit frame, it follows that x = Ql/{syh) for 0.05 < y e < 0.15! 
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Phase space cuts 




low Q 2 sample (BEMC) 


high Q 2 sample (LAr) 


Cut g 


49 GeV 2 < Q 2 < 100 GeV 2 


196 GeV 2 < Q 2 < 10000 GeV 2 


Cut g 


0.05 < y < 0.8 


Cut | 


E e > > 14 GeV 


E e , > 11 GeV 


Cut | 


157° < Oef < 173° 


30° < 9 e , < 150° 


Cut | 


20° < e q 


Cut |f 


E* > Q/10 



Table 5.2: Final phase space cuts. Note that the starred cut no. || is applied for r P , B P , 
Pe, t~c and Cp only. 



20° < e 



9qi eq. (|4.7 ), indicates the direction of the — z*-axis in the laboratory system. The 
requirement 6 q > 20° ensures a sufficient detector resolution in polar angle (s. the 
discussion in section 4.1.2). In the low Q 2 regime, this cut-off is automatically 



fulfilled owing to the previous selection in Q 2 and y\ a small number of high Q 2 



events, however, is discarded, cf. fig. 5.3 



r P , Bp, p E , t c and C P only: E* > Q/10: 

As explained in section [1]^, this cut is essential to keep the event shapes Tp, Bp, 
Pe, t~c and Cp infrared safe. In this sense, it is part of their definition. Experi- 
mentally, it ensures a minimum of hadronic activity in the CH and suppresses events 
substantially influenced by noise and leakage out of the RH. 



For an overview of the data quality cuts see table [5]3 
Data quality cuts: 

• T P , Bp, p E , t c and C P only: iV CH > 2 ; (F F) > 5-10- 5 : 

Asking for at least two objects in the CH, "unnatural" peaks at zero for pe, tq 
and Cp in the low Q 2 region are removed. "Unnatural" means that they are caused 
either by leakage out of the remnant fragmentation region or by cutting off "regular" 
jets just at the border between the two hemispheres. In fact, Bp does exhibit 
a pronounced peak at 0.5 for the rejected events showing the hadrons/clusters to 
have polar angles marginally larger than the required 90° in the Breit frame. For the 
same reason, the extreme values F min and F max are excluded by very small cut-offs. 
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£forw > 0.5 GeV: 

Requiring a minimal energy deposition in the forward region defined by 
4° < 9 < 15°, diffractive events not described by usual pQCD calculations are dis- 
carded. 

5.7° < c i < 170°: 

This cut reflects the angular acceptance for clusters that are completely contained 
in the LAr or BEMC calorimeters. 

30 GeV < (E - P z ) < 65 GeV: 

If all emerging particles of an ep reaction could be measured perfectly, then 

(E-P z ):=J2(E t - Pzt ) (5.5) 

i 

would be E p — p Zp + E e — p Ze = 2 • E e = 55 GeV. Losses due to limited acceptance, 
e.g. around the beam pipe, transform this peak into a broad distribution with tails 
down to very low values. Restricting the range of (E — P z ) effectively reduces the 
photoproduction background and the size of radiative corrections. 

low Q 2 : P t < 7.5 GeV, high Q 2 : P t < 15 GeV: 

The total transverse momentum 




of a NC event should normally be zero. To suppress remaining background from 
CC reactions or badly measured events, maximal P t 's of 7.5 GeV and 15 GeV are 
allowed in the low Q 2 and high Q 2 samples respectively. 

\(E e , - E da )/E da \ <0.25: 

According to 

E da := |g + E e (l - y da ) , (5.7) 

the energy of the scattered electron can be derived from angular information only. 
This fact is exploited e.g. for the energy calibration of the calorimeter. Asking for 
both values to be compatible within 3cr suppresses events strongly affected by 
QED radiation. 

—35 cm < z v — (z v ) < 35 cm: 

Here, it is enforced that an interaction vertex could be determined which lies within 
±35 cm ps 3cr around the average ^-position for the corresponding run period. In 
1994, (zy) ps 3 cm, and in 1995-1997, (z v ) ps -1cm. 
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Data quality cuts 




low Q 2 sample (BEMC) 


high Q 2 sample (LAr) 


Cut [7|* 


iVcH>2, F-F min >5-10" 5 , F max -F>5-1(T 5 


Cut | 


E iorw > 0.5 GeV 


Cut | 


5.7° < d < 170° 


Cut 10 


30 GeV < (E-P z ) < 65 GeV 


Cut 11 


P t < 7.5 GeV 


P t < 15 GeV 


Cut 12 


\{E e , -E da )/E da \ <0.25 


Cut tt 


—35 cm < 2 V — (z v ) < 35 cm 


Cut 14 


^backw < 10 GeV 


Cut 15 




0imp - 4>cr\ > 2° 
| ^imp ^cr | 5 CHL 



Table 5.3: Final data quality cuts. Note that the starred cut no. [7] is applied for rp, Sp, 
pp, re and Cp only. Depending on the theory to compare with, the cuts nos. ^ and |] 
may have to be considered as phase space cuts, cf. the discussion in section [7T3| . 



• ^backw < 10 GeV: 

Due to the very limited capability of the BEMC calorimeter to measure hadronic en- 
ergies, events with significant activity there are discarded. In addition, unidentified 
scattered electrons in the BEMC may lead to a rejection. 

• LAr only: |0 imp - <j> cr \ > 2°, 

-i (—65 cm < z imp < —55 cm) A -i (15 cm < z imp < 25 cm): 

In order to ensure a reliable measurement of the scattered electron, partially in- 
efficient regions such as cracks between calorimeter modules (0-cracks) or wheels 
(z-cracks) have to be avoided. imp and z imp denote the impact coordinates of the 
electron. 



The only considerable background remaining hereafter stems from photoproduction. 
In the low Q 2 sample it is estimated to be less than 3%, for the high Q 2 region it is 
negligible p3| . Residual radiative effects are accounted for by the correction procedure 
described in ch. [7|. 
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Figure 5.3: Distribution of finally selected events in the (x, Q 2 )-plane. For clarity, f 994 
data only are shown. The curves mark out the phase space cuts nos. in y, E e i, e > 
and 9 q as indicated; the cuts in y correspond to the (unlabelled) dotted lines. 



To illustrate the stability of the described selection procedure, the number of events 
accumulated in 1994 for the low Q 2 as well as the high Q 2 sample is plotted versus the 
integrated luminosity Ci nt in fig. |5]4]. Also demonstrated is the constancy of the number 
of events gathered per Ci nt of ~ lpb -1 for all four years contributing to the high Q 2 
sample. The term "pre-final" refers to the production of data n-tuples where additionally 
to the preselection the cosmic filter, the anti-photon tag and the cuts nos. §, [13] and [14] are 
in effect already. QFSELM/QESCAT is employed for the determination of the scattered 
electron. Statistics on the selection procedure can be looked up in table ^4[ 



5.5. FINAL CUT SCENARIO 
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year 


A-int/PD 


^ preselected 


^ pre-final 


jf nnai 


1994 high Q 2 


3.20 


12177 


9467 


3646 


1995 high Q 2 


4.29 


16199 


12817 


5006 


1996 high Q 2 


8.49 


30319 


23924 


9478 


1997 high Q 2 


22.26 


100466 


61711 


24477 


total high Q 2 


38.24 


159161 


107919 


42607 


1994 low Q 2 


3.18 


35051 


24690 


9761 



Table 5.4: Integrated luminosities and events gathered in the selection steps. 
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Figure 5.4: In the two plots on top, the accumulated events (upper curve: pre-final, lower 
curve: final) are drawn against Li n t for the low and high Q 2 samples of 1994. A constant 
rate of incoming events is demonstrated. The lower plot presents the constancy of the 
number of events gathered per Ci nt of ~ 1 pb _1 for all four years contributing to the 
high Q 2 sample. The statistical uncertainties are smaller than the marker sizes. 
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Chapter 6 
Event Simulation 



Measurements with the HI detector essentially comprise clusters, i.e. energy depositions in 
the calorimeters, and tracks in the tracking devices that are caused by long-lived[] particles, 
mainly hadrons. Yet, the objects dealt with in theoretical considerations of pQCD are 
partons or, equivalently, quarks, anti-quarks and gluons. Due to the complexity of the 
measuring apparatus and the underlying physical processes, a direct link from clusters 
and tracks backwards to hadrons or even partons can not be established. Nevertheless, 
information can be drawn from the selected data by comparing with model assumptions 
on a statistical basis. The first task to be carried out in this analysis chain consists in 
simulating the detector response to a given physics model. This simulation procedure is 
the subject of the next two sections. 



6.1 Simulation 
6.1.1 Parton Level 

Starting with the calculation of the matrix element of the hard scattering to 0(a 2 a s ), 
complicated integrals arise. For the purpose of generating "events," they are solved by 
employing a Monte Carlo integration technique, s. e.g. [jSSfl . Subsequently, the result has 



to be folded with pdfs describing the proton structure. For the purpose of easy access, the 
available sets of pdfs are compiled in the PDFLIB program library ||39|| . According to the 
probability distribution derived from the matrix element and the available phase space, 
a limited number of final state partons is "produced" within the framework of dedicated 
computer programs called Monte Carlo (MC) generators. For ep DIS two such programs, 
LEPTO [HOI and HERWIG El, are available. 



To account for higher orders in a Leading Logarithmic Approximation (LLA), both of- 
fer an implementation of Parton Showers (PS) including coherence effects. These branch- 
ing algorithms may be attached to incoming (Initial State PS, ISPS) and outgoing partons 
(Final State PS, FSPS) as long as their four-momentum squared t (virtuality) is above 



L Hcrc, long-lived means lifetimes t > 10 
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FSR (right) on the lepton side. 



some adjustable threshold, typically around t m 1-4 GeV 2 . Otherwise, the PS is termi- 
nated. 

Alternatively, a Colour Dipole Model (CDM) can be invoked to describe gluon ra- 
diation including the first emission in the QCDC reaction. ARIADNE |^2| supplies an 
implementation of the CDM, but is not intended to be a stand-alone program. Instead, 
it provides an interface to LEPTO so that it may be used within its framework. 

Apart from 0(a s ) QCD corrections to the Born cross section ( |3.21|) , also QED cor- 
rections of 0(a) may be sizable depending on the phase space. This is especially true 
for high y and low x. Following the singularities proportional to l/{k ■ I'), l/{k f ■ V) and 
\/{k — k' — I') 2 that appear in the calculation of the real diagrams, they can instructively 
be labelled Initial State Radiation (ISR), Final State Radiation (FSR) and Compton 
contribution, where the latter plays only a minor role. Fig. |6.1| depicts the first two of 
them where the photons are radiated predominantly collinear to the incoming respectively 
outgoing lepton. Nevertheless, all three parts are defined in the whole phase space such 
that e.g. "ISR" photons may also be directed along the scattered lepton! In addition, this 
simple picture is only valid in LLA; beyond, the separation is not unique. 

The event generator DJANGO |^3| combines the abilities of LEPTO and HERA- 
CLES j44], which provides the QED radiative effects including virtual corrections due to 
1-loop diagrams, into one software package offering the most complete description of ep 
DIS events available. 



6.1.2 Hadron Level 



Neglecting leptons and photons at this stage, we are left with partons of low virtuali- 
ties t ~ t , where perturbation theory ceases to be applicable. Lacking better theoretical 



6.2. RECONSTRUCTION 
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means, phenomenological models have to be invoked to perform the necessary fragmenta- 
tion of partons into hadrons. Since the PS cut-off to is an arbitrary parameter, the 
hadronization model should employ a similar scale from which to start in such a way that 
the dependence on i largely cancels between the PS and the fragmentation process. 

The two most popular models existing for this purpose are the string and the cluster 
fragmentation, both of which also consider colour coherence effects. The first one is 
implemented in JETSET [ 45| , which is applied for this task by the MC generators of the 
last section except for HERWIG, which makes use of the second approach. Taking into 
account decays of unstable particles as well, the final outcome of this step consists of all 
particles traversing through a real detector. 



6.1.3 Detector Simulation 

So far the involved processes were of a basically theoretical nature and could be described 
by general purpose MC programs. Since now the detector response to the passing particles 
has to be reproduced, software specifically adapted to the measuring device is required. 



The software package employed by HI is called H1SIM f46| . 

It is responsible for tracking the MC particles through a virtual HI detector and 
simulating the response signals in detail. The "events" produced that way look like real 
data. For comparison purposes it is desirable to have as many MC events as possible, at 
least about the same amount as data are available. However, the simulation step is very 
time consuming and takes several seconds of computing time per event on the computer 
systems at disposal. Therefore, the MC models to use have to be chosen carefully. 



6.2 Reconstruction 
6.2.1 Cluster Level 

Because the "raw" information (real and simulated) in the form of wire hits, cell voltages, 
etc. is not very intuitive with respect to the physics of ep scattering, the data have to 
be refined. This reconstruction process is the task of another HI software program, 
H1REC f34|, which has to be identical for real and simulated events. As a result, it 
provides i.a. particle tracks and the calorimeter clusters extensively used in this study. 
For that reason, this stage is called "cluster level." 



6.3 Comparison to Data 

Once the simulation and reconstruction have been completed, the MC models must be 
confronted with data. Ideally, they should give a good account of both, "standard" dis- 
tributions where selection cuts are applied (the energy spectrum of the scattered electron, 
dn/dE e i, for example), as well as event shape distributions dn/dF that are of special 
interest here. 
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MC model 


DJANGO 6.2 


LEPTO 6.1 


LEPTO 6.5 


HERWIG 5.8 


low Q 2 statistics 


4446 


4253 






high Q 2 statistics 


121483 


8906 


84561 


12617 



Table 6.1: Cluster level MC statistics. 



In the previous publication |23|], LEPTO 6.1 was shown to give the best description 
of the data. This is still true. However, it does not contain radiative corrections. There- 
fore, DJANGO 6.2 together with ARIADNE has been chosen as a replacement. In the 
following, the term "DJANGO 6.2" always refers to this combination! 

For testing purposes a newer LEPTO version, LEPTO 6.5, and HERWIG 5.8 have been 
used. The basic differences to the old LEPTO version are changed default settings, a new 
cut-off scheme for the divergences of the matrix elements, an improved target remnant 
treatment and the introduction of soft colour interactions to facilitate the production of 
diffractive events within DIS samples. As will be seen, this has a considerable impact on 
the event shape distributions. Statistics of the employed MC files are given in table 6.1 . In 
case of DJANGO 6.2 high Q 2 , it should be noted that an extra very high Q 2 MC file was 
produced to increase the statistics in the Q-bins seven and eight. When combined with 
the other high Q 2 data sets, it would lead to unnatural steps in some of the distributions 
presented in the next section. Therefore, it was excluded there. 



6.3.1 Standard Distributions 



In figs. [T3] and |6]4] the normalized differential distributions of E e >, 9 e >, ln(Qg/ GeV 



Uh, (E — P z ) and 9 q are presented separately for the low Q 2 and high Q 2 sample. The 
comparison with DJANGO 6.2 as well as LEPTO 6.1 reveals a good overall agreement 
with the exception of low ^-values in case of low Q 2 events. Yet, the imposed cut-off 
no. |5] of 20° < 9 q is well below the smallest occurring angle of 9 qjU1 i n ~ 30°, rendering the 
deviation harmless. 

The dips in the high Q 2 ^-distribution stem from the rejection of events with the 
scattered lepton found in areas of the LAr where the energy measurement deteriorates 
(cut no. |i~5|) . 



LEPTO 6.5 and HERWIG 5.8 are not shown here. Their agreement with data is 
similar. 



6.3.2 Event Shape Distributions 

The next step is to check the description of the event shape distributions. They are shown 



for four out of eight investigated bins in Q in figs. |3.5| and |6.6| . Here, the similarity is 



not as satisfactory as before. DJANGO 6.2 systematically tends to overshoot the data in 
the low F-region which is compensated for by underestimating them for high values in F. 
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Figure 6.2: Normalized differential distributions of Cp on cluster level. HI data (full 
symbols) are compared with LEPTO 6.5 (full lines) on the left-hand side and HERWIG 5.8 
(full lines) on the right-hand side for four out of eight investigated bins in Q. The spectra 
for (Q) = 7.5-57.9 GeV are multiplied by factors of 10 n , n = 0,1,2,3. The error bars 
represent statistical uncertainties only. 



This is worst for Tp and especially Bp, whose mean values come out too low in the MC. 
The uncorrected data means (F) are given in table p.2\ . 

Within the much lower statistics available for LEPTO 6.1, it seems to do a better job 
in reproducing the data distributions. Nevertheless, DJANGO 6.2 in combination with 
ARIADNE provides an acceptable account of the available data and, since it includes 
radiative corrections, is employed for the primary correction procedure discussed in the 
next chapter. 

As an example of a rather bad performance, fig. |6.2| presents the normalized differential 
distributions of Cp for three bins of the high Q 2 sample in comparison with LEPTO 6.5 
on the left-hand side and HERWIG 5.8 on the right-hand side. Severe deviations contrary 
to each other are observed. 
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Figure 6.3: Normalized differential distributions of E e /, 9 e / and ln(Q 2 / GeV 2 ) on cluster 
level for the low Q 2 (left) and high Q 2 (right) data sample (full symbols) in comparison 
with DJANGO 6.2 (full line) and LEPTO 6.1 (dashed line). The error bars represent 
statistical uncertainties only. 
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Figure 6.4: Normalized differential distributions of yh, (E — P z ) and 9 q on cluster level 
for the low Q 2 (left) and high Q 2 (right) data sample (full symbols) in comparison with 
DJANGO 6.2 (full line) and LEPTO 6.1 (dashed line). The error bars represent statistical 
uncertainties only. 



58 CHAPTER 6. EVENT SIMULATION 



a. 
c 

T3 



DJANGO 6.2 
LEPTO 6.1 



LU 

Q. 

C 
T3 



H1 data 

• <Q>= 7.5 GeV 
<Q) = 15.0 GeV 

a <Q> = 23.6 GeV 
<Q> = 57.9 GeV 




DJANGO 6.2 H1 data 
LEPTO 6.1 . < Q >_ 7 5 GeV 
<Q) = 15.0 GeV 
* <Q> = 23.6 GeV 
(Q> = 57.9 GeV 




0.05 0.1 0.15 0.2 0.25 

Pe 



D. 

m 

T3 



O 
T3 



10 

10 
10 
10 
10 

1 

10 " 
10 



H1 data 

• <Q>= 7.5 GeV 
<Q) = 15.0 GeV 

* (Q> = 23.6 GeV 
<Q> = 57.9 GeV 



DJANGO 6.2 
LEPTO 6.1 




0.1 0.2 0.3 0.4 0.5 



DJANGO 6.2 H1 data 
LEPTO 6.1 



• <Q>= 7.5 GeV 
<Q> = 15.0 GeV 

* <Q> = 23.6 GeV 
<Q> = 57.9 GeV 




Figure 6.5: Normalized differential distributions of the event shapes Tp, Bp, pe and tq 
on cluster level. HI data (full symbols) are compared with DJANGO 6.2 (full lines) 
and LEPTO 6.1 (dashed lines) for four out of eight investigated bins in Q. The spectra 
for (Q) = 7.5-57.9 GeV are multiplied by factors of 10™, n = 0,1,2,3. The error bars 
represent statistical uncertainties only. 
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Figure 6.6: Normalized differential distributions of the event shapes Cp, yjj and on 
cluster level. HI data (full symbols) are compared with DJANGO 6.2 (full lines) and 
LEPTO 6.1 (dashed lines) for four out of eight investigated bins in Q. The spectra 
for (Q) = 7.5-57.9 GeV are multiplied by factors of 10™, n = 0,1,2,3. The error bars 
represent statistical uncertainties only. 
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(O) 1 GeV 


(t p ) 


(Bp) 


7.46 


0.4351 ±0.0031 


0.3651 ±0.0013 


8.74 


0.4110 ±0.0034 


0.3543 ±0.0014 


14.97 


0.3456 ± 0.0021 


0.3186 ±0.0010 


17.75 


0.3203 ±0.0018 


0.3039 ± 0.0009 


23.62 


0.2782 ± 0.0020 


0.2790 ±0.0010 


36.72 


0.2367 ±0.0034 


0.2502 ±0.0019 


57.93 


0.2079 ± 0.0090 


0.2271 ±0.0053 


81.32 


0.1598 ±0.0134 


0.1891 ±0.0082 



(Q)/GeV 


(p E ) 


(tc) 


(Cp) 


7.46 


0.0797 ±0.0007 


0.1401 ±0.0013 


0.4960 ±0.0034 


8.74 


0.0805 ±0.0008 


0.1429 ±0.0014 


0.5048 ±0.0038 


14.97 


0.0737 ±0.0005 


0.1297 ±0.0008 


0.4701 ± 0.0022 


17.75 


0.0725 ±0.0004 


0.1267 ±0.0007 


0.4603 ±0.0019 


23.62 


0.0685 ±0.0004 


0.1179 ±0.0008 


0.4324 ±0.0021 


36.72 


0.0626 ±0.0008 


0.1041 ±0.0013 


0.3847 ±0.0036 


57.93 


0.0578 ±0.0021 


0.0950 ±0.0033 


0.3495 ±0.0098 


81.32 


0.0482 ±0.0030 


0.0804 ±0.0050 


0.3041 ±0.0147 



(Q) / GeV 


(Vfj) 


(Vk t ) 


7.46 


0.1273 ±0.0012 


0.1746 ±0.0023 


8.74 


0.1274 ±0.0013 


0.1471 ±0.0023 


14.97 


0.1131 ±0.0009 


0.1077 ±0.0013 


17.75 


0.1115 ±0.0008 


0.0975 ±0.0011 


23.62 


0.1047 ±0.0009 


0.0826 ±0.0012 


36.72 


0.0911 ±0.0015 


0.0603 ±0.0018 


57.93 


0.0874 ± 0.0044 


0.0483 ± 0.0040 


81.32 


0.0681 ±0.0060 


0.0341 ±0.0050 



Table 6.2: Uncorrected mean values of the event shapes and their statistical uncertainty 
as a function of Q. 



Chapter 7 

Correction Procedure 



After having established a sufficient description of data by MC models, the most important 
task remaining is disentangling the underlying physics from mere detector effects due to 
limited efficiencies and resolutions. Additionally, radiative corrections are taken into 
account. 

In principle, it is feasible to do this in a one-step procedure. Compensating effects, 
however, may mislead to the conclusion that the corrections are small, even if indeed 
they are not. Moreover, one would like to identify the dominant influence. Another 
disadvantage is the complete neglect of migrations enforced by the lack of correlations. 

For those reasons, the procedure applied here tries to differentiate between these con- 
tributions, and the unfolding of the data is correspondingly performed in two stages. The 
different correction methods available are explained in the next section. 



When only mean values do matter, the easiest thing to think of is to invoke one correction 
factor for each mean: 



7.1 Unfolding Methods 



7.1.1 Factor Method 



corr 



C-had/sim 

dat ' 




(7.1) 



7.1.2 Bin-to-bin Correction 



By applying eq. ( |7. 1| ) not to the mean values alone, but to each bin i of a differential 
distribution F i: 



1 F MC 



«, corr 




Cj,had/sim • 




(7.2) 
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the factor method can easily be extended to unfold complete distributions which again 
can be reevaluated to give a corrected mean {F) coii . 

As can be seen from eq. ( [7.2| ), however, this simple approach completely ignores the 
possibility of migrations from one bin i on hadron level to another bin j on cluster level. 
Depending on the quantity to deal with, this constitutes a severe disadvantage. A measure 
to improve the situation is to choose bin sizes in such a way that migrations are minimized. 



7.1.3 Matrix Method 



It is still better to take these migrations explicitly into consideration by employing a 



correction matrix C^. 



The matrix CV, can in principle be obtained from MC by inverting the transfer matrix T^ 
which transforms the "true" values i^ad into the observed ones 

pMC _ V^T-i pMC (7 A\ 

3 

In practice, however, the inversion leads to instabilities and oscillations unless extremely 
large MC statistics is available. Instead, we follow the strategy employed in [^3j to define 
the correction matrix to be 

Qj := (7.5) 

k 

where represents a probability density derived from MC correlations as 

j 

As a drawback of this scheme, there may remain some model dependence which can 
be overcome by an iterative procedure. 



7.1.4 Bayesian Unfolding 



The last method employs a program developed in ||47|| . It exploits Bayes' theorem on 
conditional probabilities to extract information from the observed distributions. To state 
it in a more general form, we label the fact of an event to have a "true" value F corr 
allocated to bin i as "cause" C, and the observed value Fd a t in bin j as "effect" Ej. Since 
in general the binning and also the domain of definition may be different, we have a 
number of Nc causes and Ne effects. Denoting further the initial probability for cause 
Cj to happen with Po(Cj) an d the conditional probability for the effect Ej to occur if Cj 
has already taken place with P(Ej\Ci), Bayes' theorem can be written as 

P(Cm - ^ (E>m ' P ° {Ci) . (7.7) 
E P(Ej\C k ) ■ P„(C t ) 

k=l 



7.2. CORRECTION FOR DETECTOR EFFECTS 



63 



That is, given an observation of Ej, we get a probability that it was produced by cause Cj. 

The result, however, depends on two ingredients: the conditional probabilities 
P(Ej\Ci), which can be estimated from MC, and the initial probabilities Po(Ci), that 
are a priori unknown. At first, the latter have to be assumed by the observer accord- 
ing to his prejudices which may even result in a uniform distribution in the case of 
complete ignorance. But, given a number of experimental observations with frequencies 
N(E) := {N(Ei), . . . , N(En e )}, Bayes' theorem can then be applied to gain information 
on Po(Cj). Based on these observed event numbers for the effects Ej, the expected number 
of events to be assigned to the different causes Cj can be inferred to be 

N(C t ) = ^P(C l \E 1 )-N(E 1 ). (7.8) 
i=i 

Ignoring the possibility of inefficiencies, i.e. the case when a cause does not produce one 
of the effects under consideration, the initial probability is estimated by 

P(C t ) := P(Q\N(E)) = N ^ (Cl) • (7.9) 

i=i 

If the initial distribution Po(d) is not consistent with the data, it will in general differ 
from the final one P(Cj). In an iterative procedure where the assumption on Po(Cj) is 
updated in each step, the observed events N(E) can now be exploited to unfold the "true" 
distribution N(C). 



7.2 Correction for Detector Effects 



The first step consists in deriving a correction to hadron level from the simulated 
DJANGO 6.2 MC files. Here, special care has to be taken when specifying what is meant 
by "hadron level" since there may be radiative photons. But we are interested in detector 
effects only. Thus, the main influence of these photons has to be accounted for in a later 
step. This is accomplished by first enforcing a limitation in acceptance according to cut 
no. |^, thereby throwing away radiation collinear to the incoming lepton. Furthermore, 
with regard to the finite angular resolution of the LAr calorimeter, the scattered electron 
is redefined to be merged with radiative photons if they are closer in angle than 5°. This 



value is motivated by a study in pE8fl . As a consequence, the cuts nos. pTDHT^ may reject 
events with respect to this "radiative" hadron level. 

Now one is free to draw a connection between the cluster and hadron level of the 
simulation and apply this to the data according to the four unfolding schemes described. 
Note that the selection criteria nos. |l"3|--|T5| are imposed on the cluster level only such that 
they are accounted for in this first step. 



Figs. [TTJand [T.2| show the bin-wise detector corrections as derived from DJANGO 6.2 
for four out of eight investigated bins in Q. Entries of zero (or correspondingly two, four or 
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Figure 7.1: Bin- wise detector corrections for the event shapes Tp, Bp, pe and re as 
derived from DJANGO 6.2 for four out of eight investigated bins in Q. The factors for 
(Q) = 7.5-57.9 GeV are shifted by offsets of 2 • n, n = 0, 1, 2, 3. The dashed lines delineate 
the corresponding position of unity. Entries exactly at zero (or 2, 4, 6) mean that no 
events were found on either level. The error bars represent the statistical uncertainty. 
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Figure 7.2: Bin-wise detector corrections for the event shapes Cp, yfj and yk t as derived 
from DJANGO 6.2 for four out of eight investigated bins in Q. The factors for (Q) = 7.5- 
57.9 GeV are shifted by offsets of 2 • n, n = 0,1,2,3. The dashed lines delineate the 
corresponding position of unity. Entries exactly at zero (or 2, 4, 6) mean that no events 
were found on either level. The error bars represent the statistical uncertainty. 



66 



CHAPTER 7. CORRECTION PROCEDURE 



six) occurring at the edges mean that no events were found at either level. Concentrating 
on the statistically relevant points, it must be concluded that the corrections sometimes 
are larger than agreeable. Most often they are below 30%-50%, but may rise higher in 
some bins, especially for B P and pE- 

Both, the matrix method as well as the Bayesian unfolding need the correlation ma- 
trices between hadron and cluster level as input. They are presented for each event shape 
in figs. |7.3| - [7.5| for two bins, one at low and one at high Q 2 . Note that the box sizes 
are scaled logarithmically with the number of entries. From these we can learn mainly 
two things: First, correlations, which improve with rising Q, are clearly demonstrated, 
and second, with respect to the resolution, the bin sizes are too small for the bin-to-bin 
correction. Nevertheless, the binning is kept as a compromise, because the matrix and 
the Bayes unfolding account for migrations as well as resolution effects and the corrected 
means should not be affected by too coarse bin sizes. 

With these ingredients given, two supplementary remarks are due. Concerning the 
matrix approach, the effect of iterations was revealed to be small, except for the fact that 
oscillations are gathered up. Hence, no iteration is performed here. 

With respect to the Bayesian unfolding, it is recommended in [47] to avoid building up 
large fluctuations with an increasing number of iterations by smoothing a new "initial" 
distribution P(Ci) in the intermediate steps. This proved to be difficult for the event 
shapes and finally was not done. To start from uniformly distributed initial probabilities 
did work, but required some iterations before converging, and consequently fluctuations 
showed up. It could be circumvented by choosing the MC distribution for P (Ci). Then, 
three iterations are fully sufficient. 



7.3 Radiative Corrections 

At this stage we have to perform the remaining correction from the previously defined 
radiative to a non-radiative hadron level. For that purpose a second kind of D JANGO 6.2 
MC files is necessary. They have to be generated with exactly the same settings as the 
ones before, except for the radiative effects which have to be switched off. A simulation 
is not required. Because there are no more ISR or FSR photons, the acceptance cut 
no. H may be dropped leaving the cuts nos. |10HT2l ineffective. In that way, the acceptance 
holes of the detector near the beam pipe are additionally compensated for. But due to 
our limited knowledge of that region, the extrapolation from MC possibly is not very 
precise. Notice that the cuts nos. [7] and ||| have been applied to all levels up to now and 
consequently they have not been accounted for by the unfolding procedure. Depending 
on the theory to compare with, e.g. the parton level of some MC, they may have to be 
considered as phase space cuts. Yet, they are not suitable in the framework of a pQCD 
calculation to NLO and hence are regarded as data quality cuts in this study. 

For the factor and bin-to-bin method, the combined correction steps are fully equiv- 
alent to a one-step procedure. Yet, as a matter of fact, we now have to deal with two 
different MC files. That means we do not have any correlations at our disposal such 
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that the unfolding schemes number three and four are not applicable. The bin-to-bin 
correction has to serve as the only available replacement. 

Figs. |7^ and [Z. T| display the radiative correction factors for four out of eight investi- 
gated bins in Q. Except for Tp and Bp where radiative effects are not negligible, they 
are very close to one. y^ t is special. The additionally shown histogram with the full line 
usually near "one" corresponds to the correction factors where the acceptance cut no. |] 
was kept! This demonstrates that it is not a radiative effect. yk t as applied here seems 
to be affected by the MC extrapolation, which is unfortunate and not understood at the 
moment. For all other event shapes including yfj, almost nothing changes whether cut 
no. H is active or not. In the case of Tp up to Cp this could be expected since they only 
deal with the CH opposite to the remnant in the Breit frame. 



7.4 Performance Check 

A mandatory check of correction methods is to employ one MC to unfold another. As 
a first test, this was done for models without QED radiation. LEPTO 6.1 served as the 
MC, whereas LEPTO 6.5 and HERWIG 5.8 were treated like real data. Consequently, 
the unfolded means should predominantly lie close to the hadron level means of the 
corresponding MC (almost) regardless of what LEPTO 6.1 would prefer. To be able to 
discriminate between them, the models have to make sufficiently different predictions for 
the quantity under study. From fig. |6.2| it is clear that for Cp this is certainly true. The 
outcome is presented in the upper plots of fig. [7.8| proving all four correction schemes 
to work properly. Nevertheless, all other event shapes were tested, too, yielding similar 
results, although not as evident. 

In a second test, LEPTO 6.1 events were smuggled in as data for DJANGO 6.2. Here, 
the additional complication arises that the latter does include radiative effects in contrast 
to the "data." As expected, the first correction step alone suffices since the radiative 
hadron level is designed to separate detector and radiative effects as much as possible. 
If the radiative factors are applied, too, then the sensible variables, i.e. Tp and -Bp, are 



even overdone with. This is demonstrated in the lower plots of fig. |7.8| , where the "fully" 
corrected means of the right plot deviate considerably from the hadron level means of 
LEPTO 6.1 as opposed to the plot on the left-hand side. Using only DJANGO 6.2 events 
without radiative photons, the discrepancy disappeared. 

Concluding, it can be stated that the unfolding works satisfactorily. 



7.5 Final Means 



To illustrate possible systematic effects due to the different unfolding schemes, all resulting 
corrected means are shown in figs. |7.9| and |7.10| . The statistical uncertainties encompass 
data as well as MC statistics. For the factor method, they are obtained properly by error 
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Figure 7.3: Correlation matrices for the event shapes Tp, Bp and pp as derived from 
DJANGO 6.2 for two out of eight investigated bins in Q. Note that the box sizes are 
scaled logarithmically with the number of entries. 
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Figure 7.4: Correlation matrices for the event shapes tc, Cp and y/j as derived from 
DJANGO 6.2 for two out of eight investigated bins in Q. Note that the box sizes are 
scaled logarithmically with the number of entries. 
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Figure 7.5: Correlation matrices for the event shape yk t as derived from DJANGO 6.2 for 
two out of eight investigated bins in Q. Note that the box sizes are scaled logarithmically 
with the number of entries. 



propagation. In case of an evaluation from unfolded distributions, they are estimated by 

A(F} ^ m) ivk + vk)- (7 ' 10) 

where N^at, Nuc are the appropriate event numbers and a denotes the standard deviation. 
Note that 1/ y/ A^at and 1/ V Nuc are not added conventionally in quadrature. This is 
motivated by the fact that the uncertainties derived this way are closer in value to those 
calculated for the factor method. 

For the final results the differential distributions as well as the central values are taken 
from the most sophisticated approach, i.e. the Bayesian unfolding with subsequent bin- 
to-bin radiative correction. In case of the event shape spectra, no further error evaluation 
beyond statistical ones was attempted. For a comparison with pQCD calculations, their 
presentation is postponed to ch. [H], figs. |10.1| and |10.2 . 



Concerning the mean values, the other methods serve to judge systematic effects. The 
due asymmetric uncertainty is estimated to be half of the spread caused by the maximal 
deviation to larger respectively smaller values. The sum of both then corresponds to half 
the total spread which is somewhat smaller than 2/V12 = l/y/3 as appropriate for a 
uniform distribution. 

A second source of systematic uncertainties are the electromagnetic and hadronic 



energy scales of the calorimeters. Due to recent improvements [13], the electron energies 
now can be measured to a precision of 1%, 2% and 3% up to 2i mp -coordinates of 24 cm, 
110 cm and the forward end of the LAr calorimeter respectively. Hadronic energies are 
known only to about 4%. In order to gain information on the influence of this uncertainty, 
the whole analysis is repeated scaling up and down separately for both scales. The 
discrepancies with respect to the normal central values are then attributed to two further 
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asymmetric systematic uncertainties. All three, i.e. including the unfolding, added in 
quadrature, yield the total systematics given together with the final results in table |7.1| . 

Concerning individual contributions, the impact of the electron energy, which directly 
affects the boost into the Breit frame, is almost always largest for Tp up to Cp, followed 
by unfolding effects. Realizing that hadronic systematics cancel between numerator and 
denominator for these event shapes, this is understandable. The situation is reversed for 
Ufj and yk t . Here, no cancellation occurs and the systematic uncertainty due to hadronic 
energies dominates. 



At last, figs. |7.11| and |7.12| present the final means including all uncertainties as well 
as the hadron level means of DJANGO 6.2 and LEPTO 6.1. Where available, already 
published results from [23[ and [49] are shown for comparison.^ Tp and Bp are in good 
agreement, but they should not because radiative corrections leading to smaller values 
were previously neglected. Seemingly, this was compensated by an overestimation of 
the electron energy that happened before the new calibration was introduced into the 
analysis. Since pp and r c are much less affected by QED radiation, the new data are 
shifted somewhat towards higher means. 



1 Cp is also given in 
here. 



but the values are believed to be erroneous and therefore are not reproduced 
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Figure 7.6: Bin-wise radiative corrections for the event shapes Tp, Bp, pe and tq as 
derived from DJANGO 6.2 for four out of eight investigated bins in Q. The factors for 
(Q) = 7.5-57.9 GeV are shifted by offsets of 2-n, n = 0, 1, 2, 3. The dashed lines delineate 
the corresponding position of unity. Entries exactly at zero (or 2, 4, 6) mean that no 
events were found on either level. The error bars represent the statistical uncertainty. 



7.5. FINAL MEANS 



73 



0. 

o 



10 
9 
8 
7 
6 
5 
4 
3 
2 
1 




• <Q>= 7.5 GeV 
<Q> = 15.0 GeV 

i <Q> = 23.6 GeV 
<Q> = 57.9 GeV 



. A ,.A'A--A"A-A-A~A~A~A-A-A-A'A A'A-A'A A 



_L 



0.2 0.4 0.6 



0.8 1 



10 

9 
8 
7 
6 
5 
4 
3 
2 
1 




• <Q>= 7.5 GeV 
<Q) = 15.0 GeV 

i (Q) = 23.6 GeV 
<Q> = 57.9 GeV 



> • + -f 



+ 



-A-A--A-A-'*'-A- 



a^ a -a^a^a«ai a i; a ; 



s 


10 


>; 




T3 


9 


ra 
.c 




■3 
ra 


8 


^ 




c 

o 


7 




6 




5 




4 




3 




2 




1 








• <Q>= 7.5 GeV 
<Q) = 15.0 GeV 

i (Q) = 23.6 GeV 
<Q> = 57.9 GeV 



0.2 



0.4 



0.6 




Figure 7.7: Bin-wise radiative corrections for the event shapes Cp, yfj and y\. t as derived 
from DJANGO 6.2 for four out of eight investigated bins in Q. The factors for (Q) = 7.5- 
57.9 GeV are shifted by offsets of 2 • n, n — 0, 1, 2, 3. The dashed lines delineate the 
corresponding position of unity. Entries exactly at zero (or 2, 4, 6) mean that no events 
were found on either level. The error bars represent the statistical uncertainty For 
yk t the additional histogram (full line) shows the correction factor without acceptance 
extrapolation. 
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Figure 7.8: Top: Unfolded means of Cp as a function of Q for simulated high Q 2 
LEPTO 6.5 (left) and HERWIG 5.8 (right) event samples. LEPTO 6.1 was employed 
as unfolding MC. Bottom: Unfolded means of Tp as a function of Q for simulated 
LEPTO 6.1 event samples with (right) and without (left) radiative corrections derived 
from DJANGO 6.2. In addition, the corresponding hadron levels of the MC's are shown. 
The error bars represent statistical uncertainties only. 
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Table 7.1: Corrected mean values of the event shapes as a function of Q. The first 
error is statistical, the second systematic and comprises unfolding as well as energy scale 
uncertainties added quadratically. 
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Figure 7.9: Unfolded means of Tp, Bp, pe and tc as a function of Q for the four described 
correction schemes. DJANGO 6.2 served as unfolding MC. Note that the step from 
radiative to non-radiative hadron level can be accomplished by the factor or bin-to-bin 
methods only. The latter was employed for this purpose in case of the matrix and Bayes 
procedures. The error bars represent statistical uncertainties only. 
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Figure 7.10: Unfolded means of Cp, t/fj and yk t as a function of Q for the four described 
correction schemes. DJANGO 6.2 served as unfolding MC. Note that the step from 
radiative to non-radiative hadron level can be accomplished by the factor or bin-to-bin 
methods only. The latter was employed for this purpose in case of the matrix and Bayes 
procedures. The error bars represent statistical uncertainties only. 
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Figure 7.11: Corrected means (full symbols) of Tp, Bp, and tq as a function of Q. 
The inner error bars represent the statistical uncertainty, the outer ones the total statistic 
and systematic uncertainty. They are compared with the hadron level mean values of the 
DJANGO 6.2 (dashed) and LEPTO 6.1 (dotted) MC and with already published results 
by HI (hollow symbols). 
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Figure 7.12: Corrected means (full symbols) of Cp, yjj and y^ t as a function of Q. The 
inner error bars represent the statistical uncertainty, the outer ones the total statistic 
and systematic uncertainty. They are compared with the hadron level mean values of the 
DJANGO 6.2 (dashed) and LEPTO 6.1 (dotted) MC. 
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Chapter 8 

NLO Integration Programs 



According to the definition of the event shapes F in section [4.2| , they will take on non- 
trivial values, i.e. F ^ 0, when derived to 0(a s ) at least. In terms of cross sections, the 
distributions dn/dF may be written to LO accuracy as 

' C 1 (F)a s . (8.1) 



0tot dF 



The first moment or mean value, characterizing the distribution in F of an event sample 
in a simple way, now reads 



■ max 



1 F ^ dF 1 F T da 
(F):= -J = ± / FpF, (8.2) 

-fmax (J. . I 

I %dF { 
o 

which can be expanded to LO in a s as well: 

(F) = c liF a s . (8.3) 

However, theoretical arguments show that the achieved precision of LO calculations is 
not sufficient to extract much information on pQCD from experimental data. Especially 
the concept of a running coupling constant a s (/J, a ), explained in the next section, requires 
the distributions and mean values to be known at least to NLO, i.e. 0(a 2 s ). 



8.1 The Running Coupling Constant 

Following an approach from ||50|| , one may assume that a dimensionless quantity (F) de- 



pends on a single energy scale Q which is much larger than any other parameter measured 
in units of energy such that quark masses may be set to zero.Q Then, the conclusion is 

1 Note that this step requires a more careful consideration, cf. the ref. given above. 
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that (F) has a constant value independent of Q; there is no other possibility to remove 
the dimension of Q. However, this is not true in a renormalizable quantum field theory 
such as QCD. 

Perturbative calculations operating with "naked" charges eo (or masses mo) of point 
particles have to deal with ultraviolet divergences occurring in the integrations of loop 
momenta d A k /(27r) 4 . The procedure allowing us to absorb these divergences in physical 
charges e is called renormalization. Results concerning the (measured) physical charges 
are then all finite. 

There exist several recipes how this can be done [|51 . One of the most popular, 
also adopted in this analysis, is the modified minimal subtraction scheme MS, where the 
singularities are regulated by a reduction of the space-time dimensions to n < 4: 



d A k 



fi 



2c 



d A ~ 2e k 



(2*) 



4-2e 



U) 



At an arbitrary scale called renormalization scale fi R , the integrals are subdivided into 
finite and divergent parts for e — > 0. The singular terms (and some finite ones depending 
on the renormalization scheme) are then attributed to a renormalized charge e. However, 
the remaining finite terms also modify the charge e or equivalently the coupling constant 
a := e 2 /47r, which now explicitly varies with fi R \ 

Relating this to QCD, we have a strong coupling constant a s := g 2 /^ containing 
fi R - - a second mass scale that enters the game. (F) may therefore depend on the 
ratio Q 2 / fi 2 R . But since fi R is an arbitrary parameter of the renormalization procedure, it 
should not make an appearance in physical observables like (F). This can be expressed 
mathematically by 



or, using the notation 



2 d 2 & a s 9 



In 



/4 



(5(a s 



Ma 



2 da s 



d/4 



(F)(Q 2 /ft 2 R ,a s )=0 (8.5) 

(8.6) 
(8.7) 



as 



9 n , \ 9 

-dt+^dcV. 



(F) (e\a s 







18) 



A solution to this partial differential equation can be given implicitly by defining the 
running coupling a s {Q 2 ) as 



dx 

W) 



a. 



19) 
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such that (F) (l,a s (Q 2 )) fulfils eq. QB.Bp. The scale dependence of any (F) calculated 
perturbatively now follows from a s (Q 2 ) provided eq. Q8.9|) could be solved. Since QCD is 
an asymptotically free theory, i.e. lim a s (Q 2 ) = 0, this can also be achieved by expanding 

Q— >oo 

the /3-function of the Renormalization Group Equation ( |8.7| ) (RGE) in powers of a s : 

n=0 

The first two coefficients defined that way are 

ft = (8.11) 

ft = ,8,2) 

where Nj is the number of active flavours. /?2 and further coefficients are scheme depen- 
dent. 

Retaining only the first term, the 1-loop solution for a s (Q) can be written as 

<Q) = Z asiflR) ,„^ ■ (8-13) 



1 + ga.W In (£) 



Using this equation, a s can be evaluated at any (sufficiently high) scale provided it is 
known at one point, say Mz- Historically, a dimensional parameter A derived from 

oo 

ln x^ = - / ( 8 - 14 ) 

A 2 J (3{x) 

was determined instead of a s (Mz). The 1-loop formula then reads 

as{Q) = TTJo\ ( 8 - 15 ) 

and demonstrates that A indicates the scale at which the coupling would diverge if ex- 
trapolated outside the perturbative regime. 

In principle, both approaches are equivalent and any measurement of a s (Mz) can 
be converted into a corresponding value for A and vice versa. The latter, however, has 
some disadvantages: It is not dimensionless, it depends on Nf and on the renormalization 
scheme^ and — most important — there exist two slightly different definitions for A in 
the literature || ^0 . 



2 That is, it has to be labelled e.g. A 



ATOMS' 
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Except for the consistency check in section |9.2| employing A 5 jjjg, we therefore follow 
ref. |52| and apply the 2-loop equation according to 

fris a s {Mz) 

a s {Q) = 77T7T\ ( 8 - 16 ) 

1 + a s (M z ) ■ (^) 



where 



L {1) (^-) =-ln#- (8.17) 
1 MzJ 2tt M z 



reproduces the 1-loop formula above and 



*(£)-(£+£w*>>£ 



gives the 2-loop result. 

Returning to (F), the perturbative expression to all orders 



(F) = X>W<W (8-19) 

n=l 

does not depend on \i R . But we do not know all coefficients, and therefore any approxi- 
mation by a truncated series does. To be more precise, 

d N 

^TT E ^J«?W « +1 (^)) • (8-20) 

"^fl n= l 

According to current knowledge, the dependence on fi R is much reduced if more terms are 
inserted. Hence, the theoretical uncertainty usually associated with the choice of scale, 
whose variation reflects the unknown higher orders, is diminished. Up to now, we merely 
considered ci, which itself does not contain fi R , and the complete /^-dependence stems 
from a s (fi R ) and is rather large. It is reduced with the inclusion of C2{ijl r )ol 2 s {ijl r ) with 

c 2 (Q) = c 2 (l) + § In ^ci, (8.21) 

which can be calculated by the three NLO integration programs presented in the next 
sections. 



8.2 NLO Integration Techniques 

As mentioned in section |3.4| on the cross section to 0(a s ), again infrared and collinear 
divergences are present in calculations to 0(a 2 s ) which cancel between real and virtual 
diagrams. Examples of real corrections to the QCDC and BGF processes with now three 
final state partons are given in fig. |8.1| . 
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Figure 8.1: Examples of real corrections to QCDC (left) and BGF (right) Feynman graphs. 



In order to eliminate the occurring singularities and obtain numerically stable results, 
the necessary cancellations were performed analytically for a certain jet algorithm in the 
programs DISJET [53| and PRO JET J54|. However, these programs lack the flexibility to 
evaluate NLO corrections for arbitrary infrared-safe observables and they contain approxi- 
mations that are inappropriate for certain phase space regions. Therefore, more elaborate 
methods have been developed that ensure the cancellations and evaluate numerical results 
for arbitrary observables. Basically, two such schemes are known: the phase space slicing 
method |155| and the subtraction method 



In the phase space slicing method, the separation of singular from non-singular phase 
space regions is done by surrounding partons with a small cone. Real emissions of partons 
outside the cone defined by = (pi + Pj) 2 > s m m, where s m i n ~ 0.1 GeV 2 , are computed 
exactly. Inside the cone, soft and collinear approximations are used. A drawback of this 
procedure is the possibility of a residual dependence on the technical parameter s m i n . For 
any calculation it has to be established that a plateau has been reached for small enough 
cut-offs. In addition, this scheme prevents the integration of the event shape distributions 
Fda/dF to be carried out over the whole phase space from zero to -F max . Instead, a lower 
bound F cut > has to be imposed and s min has to be adjusted accordingly. Moreover, the 
introduction of crossing functions |5"E |, that have to be evaluated for every parton density 
parameterization, is required. 

In case of the subtraction method, the singularities are cancelled point-by-point in 
phase space. A cut-off like s min is not formally required by this scheme, although a tiny 
value of s min ~ 10~ 8 -10~ 10 GeV 2 is applied to set the number of significant digits. Hence, 



the event shape distributions Fda/dF can be integrated down to zero [p7 



The first method is employed in the NLO integration program MEPJET H58H , whereas 
the subtraction procedure is used in DISENT |59[ and DISASTER++ pD[ . All three 
programs, however, do not contain the full set of diagrams to 0(a 2 ). Additional 2- 
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Pi 



^p 



pi 



i/p 



M/p 



Figure 8.2: Examples of 2-loop virtual correction Feynman graphs. 



loop virtual corrections from graphs like the ones in fig. |8]^ appear and are difficult to 
handle. Yet, they contribute as interference term to the Born process only and therefore 
do not affect the event shape distributions except for the bin at zero. The integral 
J max FjpdF still contains integrable divergences, but they are not harmful for DISENT 

or DISASTER++ |TT|. A calculation of the total cross section J Fmax jpdF, needed for 
the normalization, to 0(a 2 ) (i.e. NNLO), however, is not possible. A Taylor expansion of 
eq. (pi) 



(F) 



r max 

o 



—dF 

dF ar 



r max 

/ 

o 



da dF 



aia\ + a 2 a 2 
b a Q + biot, 1 + b 2 a 2 



ai ! 



a 2 - 



qih 
bo 



a 2 + 0(al 



(8.22) 

(8.23) 
(8.24) 



A) 



shows that, owing to the property of F to be zero in the QPM limit and hence a® 
Otot needs only be known to 0(a s ) for (F) to be still correct to 0(a 2 ). 



0. 



8.3 First Comparisons of NLO Programs 

In view of the complexity of the calculations and their implementations into program code, 
it is mandatory to compare the results where possible. A first comparison of MEPJET 1.4 
and DISENT 0.0 was performed in [p5fl . Due to the phase space slicing method applied 
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Figure 8.3: Differential distributions of Tp for four out of seven investigated bins in Q. 
LO (left) and NLO (right) pQCD calculations of MEPJET 1.4 (hollow symbols) and 
DISENT 0.0 (full lines) are compared. The spectra for (Q) = 7.5-62.7 GeV are multiplied 
by factors of 10 n , n = 0, 1, 2, 3. 



in MEPJET, the differential distributions for F > F cut > only could be considered. 
In order to keep F cut small, a cut-off of s min = 0.01 GeV 2 instead of 0.1 GeV 2 was used. 
Nevertheless, the left-most bins extending down to zero should not be taken seriously, even 
if they accidentally do agree. Fig. |873| shows a perfect agreement of the LO spectra and, 
within the calculational precision of 3-5%, compatible results for the NLO distributions. 
With the advent of DISASTER++, a check of all three programs against each other 



was carried out and some discrepancies have been revealed [p0(1 . Partially, these are 
understood and meanwhile an improved version DISENT 0.1 has been released. In event 
shape spectra like fig. |S.3| , effects are hardly visible, but the mean values have increased 
considerably in the low Q 2 region compared to the calculations that entered in The 
worst case was the jet broadening Bp, whose mean value at Q ~ 7.5 GeV rose by about 
14%. 

The most recent program versions are MEPJET 2.2, DISENT 0.1 and DISAS- 
TER++ 1.0.1. In this analysis only DISENT will be used, but more elaborate com- 
parisons are under way in the DESY workshop on Monte Carlo Generators for HERA 
Physics. 
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Chapter 9 
DISENT Results 



DISENT, the NLO integration program mainly employed here, is structured similarly 
to a MC event generator. A main routine generates configurations of parton momenta 
in the allowed phase space region — for our purposes defined by the cuts nos. of 
section 5J5 — , stores them in an event record and calculates the weight according to the 
matrix elements. Consequently, a user routine may analyze these events with respect to 
arbitrary infrared and collinear safe observables. For physical cross sections, the weights 
have to be multiplied by the strong coupling constant raised to the proper power and 
by weights extracted from the pdfs of the proton. These are not inherent to DISENT 
but must be gathered from somewhere else, e.g. the PDFLIB [3E]. As standard options 
we use the set MRSA'-115 and Nf = 5 as number of active flavours. Note that 
DISENT currently neglects quark masses as well as electroweak corrections proportional 
to Q 2 1 \Q 2 + Mf) and (Q 2 /(Q 2 + Mf)) 2 . Fortunately, only the highest Q-bin investigated 
with (Q) = 81.32 GeV may be affected noticeably |I3] |. In the near future, DISENT may 
be extended to include the electroweak terms |63 |. 

In distinction to MC generators, however, several different but correlated event records, 
in total 14, are produced for one DISENT "event." They correspond to contributions 
of real and virtual diagrams plus counter terms that ensure the necessary cancellation 
of divergences. These contributions have coupled weights that are positive as well as 
negative. For the error evaluation, it is mandatory that these weights are allowed to 
counterbalance each other before finally e.g. being histogrammed. This may be done by 
creating an intermediate array or histogram, where only the contributions of one event 
are entered. All bins not empty at the end are then transferred to the final array with 



their respective weights. A more detailed program description can be found in [54 



9.1 Evaluation of the Statistical Uncertainty 

Given the framework of the DISENT program, the calculation of differential distributions 
da/dF and mean values (F) is straightforward. Positive and negative weights of the 
different contributions are summed up. Yet, the evaluation of the statistical uncertainties, 
where the technique of intermediate histograms is presupposed, is more involved. 
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First, we want to consider cross sections that can be written as a sum of weights Wi, 
given by DISENT, which are already normalized to the number of events N: 



N 



a = Wi . 



(9.i; 



i=i 



By reintroducing unnormalized weights Vi, one obtains 



a 



1 N 

8=1 



(9.2) 



demonstrating that a is computed analogously to a mean value. Its statistical uncertainty 
can be estimated by 



- 1 I iV^ V< A^ 2 I ^ 

\ t=l \i=l 



(9.3) 



\ 



N 



y i=l \i=l 



With iV « JV — 1, this finally gives 



A 



(7 



iV 

2 V=l 



AT 



(9.4) 



which in the limit of very large N becomes 



Aa 



N 



For our purposes the second term is not negligible and formula (]9. 
Turning to the event shape means 

N 

F i W i 

i=l 



( F ) v 



(9.5) 

has to be applied. 

(9.6) 



i=i 
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one could be misled to the assumption that a similar uncertainty as in eq. ( |9.4p could be 
derived here. Yet, formula ( |9.6|) emphasizes that the weights Wi and w i are not identical 
because the numerator is evaluated to 0(a 2 s ) and the cross section in the denominator to 
0(a s ) only. Instead, the uncertainties have to be estimated separately, and afterwards 
the error of the quotient can be obtained by the usual formulae. The final result reads: 



A (F) = (F) 



\ 



N N 

E ( W E < 2 

i=l , i=l 



N 

E F i W > 
i=l 



N 



E w i 



N 



(9.7) 



As a cross-check, the same DISENT 0.0 calculation has been repeated about a hundred 
times with different seeds of the random number generator for 0.1, 0.2, 0.5 and 1.0 million 
generated parton events each. Fig. |9.1| shows on the left-hand side histograms of the 



obtained cross sections cr to t and mean values (rp) together with the fit of a Gaussian dis- 
tribution for half a million events. It describes the distribution of the results well and its 
standard deviation is compatible with the deduced uncertainty RMS ca i c . The right-hand 
side presents the dependence of the uncertainties for both, the statistical and the compu- 
tational derivation, on the number of events produced. It nicely exhibits errors shrinking 
proportional to l/viV. The computation underestimates only slightly the statistically 
achieved variance. 

The error bars in differential distributions are evaluated according to eq. 



9.2 x -Dependence 

a s {Mz) is the fundamental parameter of QCD that has to be determined by experiment. 
An assumption on it, however, is already contained in our cross section and mean value 
calculations. Therefore, they cannot directly be used to extract information from our 
data. In |25| and |23| this problem was avoided by fitting eq. ( p.l9j ) up to 0(a 2 s ) to 



the obtained mean values of each event shape. The resulting coefficients c%,f an d C2,f 
are subsequently employed for an a s (M^)-dependent parameterization of the event shape 
means. 

In contrast to e + e~ physics, where c^p and C2,f truly are pure numbers (s. e.g. [p5]), 



this is not the case in ep DIS because in addition to the matrix elements x-dependent 
pdfs enter their evaluation. In general, the coefficients will therefore be functions of x 



which was taken into account in an approximate way only by the above ansatz. Fig. £)72 
presents as examples for a strong — the "worst" case being y^ t — and a weak variation 
of the mean values with x curves of (rp) and (Cp) versus Q for four different bins in x. 
Note that in this investigation the two highest Q-bins were merged. 

A solution to this problem is the direct calculation of c\ t F and C2,f for every Q-bin 
with its specific x-range separately. This method corresponds to an approximation of 
functions c nt p(x) by step functions instead of one suitably chosen number. 



92 



CHAPTER 9. DISENT RESULTS 




Figure 9.1: Left: Histograms of <T to t (top) and (rp) (bottom) plus Gaussian fit resulting 
from 117 DISENT 0.0 calculations with half a million events each. RMS ca i c denotes the on 
average estimated standard deviation. Right: Dependence on 1/y/N of the statistical (full 
symbols) and computational (hollow symbols) uncertainty estimates for cr to t and (rp). 




Figure 9.2: Mean values of rp (left) and Cp (right) versus Q in four different bins of x. 



9.3. FINAL PERTURBATIVE COEFFICIENTS 



93 




Figure 9.3: Fits (dashed and full lines) of A 5 j^g for the LO (hollow symbols) and NLO 
(full symbols) means of rp (left) and Cp (right). The steps express the x-dependence of 
the perturbative coefficients c^p and C2,f- 

As a consistency check, the perturbative mean values at LO and NLO are fit with the 
coefficients derived in the same computation. The outcome for A 5 should be identical 
to the initial choice, that is A 5 jjjg = 179 MeV for MRSA'-115. The achieved fits for (rp) 
and (Cp) are displayed in fig. |9~3| . Again, the much steeper steps for (t p ) indicate its 
stronger x-dependence. The A 5 ^§ values extracted turn out to be very well compatible 
with the expectation for all event shapes except for the jet broadening Bp in NLO! This 
may be a hint that there is a problem in the DISENT 0.1 calculations at low Q 2 . 



9.3 Final Perturbative Coefficients 



The final results obtained with DISENT 0.1 are compiled in the tables |9.1| - |9.8| . The first 
one contains the mean Q- and x-values to NLO for the eight bins of section |0| as well as 
the LO and NLO cross sections. The latter are differentiated by indices "1H" and "2H" to 
indicate whether the CH alone (1H, energy cut no. |6]) or both hemispheres (2H, no energy 
cut) are allowed to contribute to the hadronic final state. To LO there is no difference 
because the only parton available has E* = Q/2 and runs into the CH by definition. 

Tables p^2]-p78| comprise the mean values of all event shapes for our standard pdf 
set MRSA'-115 and the perturbative coefficients c^p and C2,p needed for the a s (In- 
dependent parameterization of the means. 
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(Q) 1 GeV 


(x) 


crLo/pb 




0"nlo,2h/p& 


7.46 


4.58 • 10~ 3 


2519.9 ±0.3 


2232.5 ±0.4 


2402.6 ± 0.4 


8.80 


7.21 • 10~ 3 


2357.5 ± 0.4 


2111.9 ±0.5 


2245.8 ±0.5 


14.95 


1.60 • 10~ 2 


534.4 ±0.1 


484.4 ±0.1 


503.8 ±0.1 


17.73 


2.16- 10~ 2 


584.7 ±0.1 


533.4 ±0.1 


551.1 ±0.1 


23.75 


3.16- 10- 2 


465.8 ±0.1 


428.4 ±0.1 


438.9 ±0.1 


36.69 


5.80- 10~ 2 


156.2 ±0.1 


145.9 ±0.1 


147.9 ±0.1 


57.61 


1.09- 10" 1 


24.79 ± 0.02 


23.63 ±0.02 


23.75 ±0.02 


80.76 


1.74- 1Q- 1 


7.502 ±0.006 


7.291 ±0.006 


7.305 ± 0.006 



Table 9.1: Mean Q- and x-values as well as the LO and NLO cross sections as calculated 
by DISENT 0.1. The indices "1H" and "2H" differentiate between cross sections with 
and without the minimal energy cut in the CH. 



(Q) / GeV 


( t p)lo 


( t p)nlo 


c l,rp 


C 2,r P 


7.46 


0.2633 ±0.0001 


0.3377 ±0.0014 


1.4559 ±0.0006 


2.250 ±0.041 


8.80 


0.2391 ±0.0001 


0.3070 ±0.0012 


1.3735 ±0.0008 


2.226 ±0.038 


14.95 


0.1847 ±0.0001 


0.2387 ±0.0009 


1.1895 ±0.0008 


2.222 ±0.038 


17.73 


0.1684 ±0.0001 


0.2173 ±0.0008 


1.1217 ±0.0007 


2.159 ±0.037 


23.75 


0.1453 ±0.0002 


0.1871 ±0.0013 


1.0209 ±0.0011 


2.065 ±0.062 


36.69 


0.1138 ±0.0001 


0.1480 ±0.0007 


0.8622 ±0.0011 


1.978 ±0.041 


57.61 


0.0846 ± 0.0002 


0.1113 ±0.0008 


0.6909 ±0.0014 


1.811 ±0.055 


80.76 


0.0654 ±0.0002 


0.0864 ±0.0006 


0.5622 ±0.0022 


1.604 ±0.043 



Table 9.2: Mean values and perturbative coefficients for t p as derived from DISENT 0.1. 



(Q) / GeV 


( B p)lo 


( B p) NLO 


C 1,B P 


C 2,B P 


7.46 


0.3923 ±0.0001 


0.2524 ±0.0018 


2.1696 ±0.0007 


-3.444 ± 0.049 


8.80 


0.3604 ±0.0002 


0.2510 ±0.0018 


2.0699 ±0.0009 


-2.894 ±0.052 


14.95 


0.2879 ±0.0001 


0.2309 ±0.0014 


1.8545 ±0.0009 


-1.847 ±0.051 


17.73 


0.2659 ±0.0001 


0.2215 ±0.0012 


1.7713 ±0.0008 


-1.513 ±0.047 


23.75 


0.2350 ±0.0002 


0.2063 ±0.0013 


1.6515 ±0.0013 


-1.034 ±0.056 


36.69 


0.1921 ±0.0002 


0.1795 ±0.0009 


1.4565 ±0.0015 


-0.437 ±0.047 


57.61 


0.1521 ±0.0002 


0.1528 ±0.0007 


1.2431 ±0.0021 


0.245 ± 0.042 


80.76 


0.1254 ±0.0002 


0.1309 ±0.0005 


1.0798 ±0.0023 


0.571 ±0.034 



Table 9.3: Mean values and perturbative coefficients for Bp as derived from DISENT 0.1. 
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(Q) / GeV 


(Pe)lo 


(Ps)nlo 






7 46 


D D464 + 0001 

U.UtiUt: _1_ W.WWWX 


0538 ± 0006 


2564 ± 0002 


D 243 ±0017 

V . Z- t: - J _1_ W.WX 1 


8.80 


0437 ± 0001 


0522 ± 0005 


2509 ± 0002 


293 ± 014 


14.95 


0.0367 ±0.0001 


0.0464 ±0.0004 


0.2363 ±0.0002 


0.405 ±0.015 


17.73 


0.0347 ±0.0001 


0.0445 ± 0.0003 


0.2309 ±0.0003 


0.434 ±0.014 


23.75 


0.0316 ±0.0001 


0.0406 ±0.0009 


0.2222 ±0.0004 


0.448 ± 0.040 


36.69 


0.0274 ±0.0001 


0.0367 ±0.0003 


0.2076 ±0.0004 


0.537 ±0.019 


57.61 


0.0233 ±0.0001 


0.0321 ±0.0004 


0.1907 ±0.0007 


0.589 ±0.028 


80.76 


0.0204 ±0.0001 


0.0279 ±0.0002 


0.1761 ±0.0009 


0.569 ±0.017 



Table 9.4: Mean values and perturbative coefficients for pe as derived from DISENT 0.1. 



(Q) / GeV 


( t c)lo 


( t c)nlo 


C l,T C 


C 2 ,r c 


7.46 


0.0701 ±0.0001 


0.0907 ±0.0010 


0.3876 ±0.0003 


0.620 ±0.028 


8.80 


0.0664 ±0.0001 


0.0872 ±0.0008 


0.3816 ±0.0004 


0.675 ±0.025 


14.95 


0.0569 ±0.0001 


0.0770 ±0.0007 


0.3666 ±0.0004 


0.814 ±0.026 


17.73 


0.0542 ±0.0001 


0.0738 ±0.0006 


0.3607 ±0.0004 


0.854 ±0.025 


23.75 


0.0500 ±0.0001 


0.0688 ±0.0006 


0.3515 ±0.0006 


0.911 ±0.030 


36.69 


0.0443 ± 0.0001 


0.0618 ±0.0005 


0.3357 ±0.0006 


0.999 ±0.030 


57.61 


0.0389 ±0.0001 


0.0550 ±0.0005 


0.3179 ±0.0007 


1.077 ±0.034 


80.76 


0.0348 ±0.0001 


0.0490 ±0.0004 


0.2991 ±0.0009 


1.073 ±0.028 



Table 9.5: Mean values and perturbative coefficients for r c as derived from DISENT 0.1. 



(Q) / GeV 




(C^NLO 


c i,c P 


°2,C P 


7.46 


0.2708 ± 0.0002 


0.3470 ±0.0035 


1.4977 ±0.0010 


2.304 ±0.096 


8.80 


0.2570 ± 0.0002 


0.3325 ±0.0031 


1.4765 ±0.0013 


2.472 ±0.091 


14.95 


0.2207 ±0.0002 


0.2945 ±0.0024 


1.4220 ±0.0015 


3.000 ±0.094 


17.73 


0.2103 ±0.0002 


0.2826 ±0.0021 


1.4008 ±0.0015 


3.151 ±0.088 


23.75 


0.1945 ±0.0003 


0.2638 ±0.0024 


1.3679 ±0.0022 


3.371 ±0.112 


36.69 


0.1728 ±0.0003 


0.2374 ±0.0019 


1.3109 ±0.0024 


3.689 ±0.107 


57.61 


0.1522 ±0.0004 


0.2126 ±0.0019 


1.2451 ±0.0033 


4.049 ±0.124 


80.76 


0.1371 ±0.0004 


0.1909 ±0.0013 


1.1796 ±0.0044 


4.063 ±0.095 



Table 9.6: Mean values and perturbative coefficients for Cp as derived from DISENT 0.1. 
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(Q) / GeV 


(y/j) L o 


(Vfj)mx> 






7.46 


0.1694 ±0.0001 


0.1575 ±0.0007 


0.9369 ±0.0004 


-0.197 ±0.020 


8.80 


0.1553 ±0.0001 


0.1493 ±0.0007 


0.8918 ±0.0005 


-0.043 ±0.021 


14.95 


0.1219 ±0.0001 


0.1253 ±0.0006 


0.7852 ±0.0005 


0.259 ±0.021 


17.73 


0.1119 ±0.0001 


0.1172 ±0.0005 


0.7452 ± 0.0005 


0.341 ±0.020 


23.75 


0.0976 ±0.0001 


0.1054 ±0.0005 


0.6859 ±0.0007 


0.473 ± 0.024 


36.69 


0.0779 ±0.0001 


0.0875 ± 0.0004 


0.5905 ±0.0007 


0.618 ±0.019 


57.61 


0.0594 ±0.0001 


0.0697 ±0.0003 


0.4854 ±0.0007 


0.739 ±0.018 


80.76 


0.0470 ±0.0001 


0.0564 ±0.0002 


0.4052 ±0.0009 


0.739 ±0.014 



Table 9.7: Mean values and perturbative coefficients for yfj as derived from DISENT 0.1. 



(Q) / GeV 


(Vkt)u> 


(Vkt) NLO 


C l,Vk t 




7.46 


0.2138 ±0.0001 


0.2942 ±0.0004 


1.1825 ±0.0004 


2.529 ±0.016 


8.80 


0.1803 ±0.0001 


0.2366 ±0.0003 


1.0348 ±0.0005 


1.929 ±0.016 


14.95 


0.1214 ±0.0001 


0.1514 ±0.0002 


0.7818 ±0.0004 


1.299 ±0.013 


17.73 


0.1045 ±0.0001 


0.1288 ±0.0002 


0.6960 ±0.0004 


1.123 ±0.011 


23.75 


0.0830 ±0.0001 


0.1008 ±0.0002 


0.5830 ±0.0005 


0.920 ±0.013 


36.69 


0.0562 ±0.0001 


0.0676 ±0.0001 


0.4256 ± 0.0004 


0.686 ±0.009 


57.61 


0.0349 ±0.0001 


0.0423 ±0.0001 


0.2851 ±0.0007 


0.515 ±0.008 


80.76 


0.0235 ±0.0001 


0.0287 ±0.0001 


0.2015 ±0.0009 


0.405 ± 0.005 



Table 9.8: Mean values and perturbative coefficients for y^ t as derived from DISENT 0.1. 



Chapter 10 
Power Corrections 



To get a first impression of the extent to which pQCD is able to approximate the data, 
figs. 10.1 and |10.2| present the normalized differential distributions of the seven inves- 
tigated event shapes, corrected to hadron level, in comparison with DISENT 0.1 NLO 
calculations. Note that the first data bin was left out in both normalizations. This is 
necessary since according to section |8]2] it is not possible to go down to F = in 0(a 2 s ). 
Starting with bin one, the NLO histogram could alternatively be normalized to ctnlo 
yielding a slight mismatch by a factor of ct NL o/o"nnlo- Except for a small shift in height 
complicating a comparison of shapes, the conclusions would not change. 

Obviously, there are large discrepancies between data and the NLO predictions. Com- 
mon to all event shapes, however, is the clear tendency that the description of data 
improves considerably with rising Q; to what degree depends on the variable under con- 
sideration, yfj and yk t exhibit a very fast convergence hinting at small hadronization 
corrections, whereas the prediction for Cp approaches the data slowly. 

Turning back to the mean values (F), we can, except for the y shapes, expect them 
to be substantially larger than predicted by pQCD with decreasing differences for in- 
creasing Q. Indeed, physical quantities receive, in addition to pQCD, non-perturbative 
contributions which basically behave power-like, i.e. oc 1/Q P , hence called power correc- 
tions. The cross section e + e~ — » hadrons, o~h.ad-> for example exhibits a l/Q 4 -term which 
is also understandable on theoretical grounds |j(J. Less inclusive variables like event 
shapes are not as well-behaved as <7h a d and reveal 1 / Q-corrections that are sizable even at 
Q — Mz They can be made plausible in the framework of a longitudinal phase space 



or tube model ||68|| , which essentially is the simplest version of the string fragmentation im- 
plemented in JETSET. One possibility to compare with pQCD is based on hadronization 
models like the one implemented in JETSET by unfolding the data in a third correction 
step to some kind of (ambiguous) parton level. 

We do not follow this recipe. Instead, we will call upon explicit formulae trying to 
parameterize these effects. They are presented in sections |10.1| and |10.3| for the tube model 
and a new approach respectively. The main goal of this work is to test these models. 
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Figure 10.1: Normalized differential distributions of the event shapes r P , B P , p E and 
Tc corrected to hadron level. HI data (full symbols) are compared with DISENT 0.1 
NLO calculations (full lines) for four out of eight investigated bins in Q. The spectra 
for (Q) = 7.5-57.9 GeV are multiplied by factors of 10™, n = 0,1,2,3. The error bars 
represent statistical uncertainties only. Note that the first data bin was left out in both 
normalizations. 
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Figure 10.2: Normalized differential distributions of the event shapes Cp, yjj and y^ t 
corrected to hadron level. HI data (full symbols) are compared with DISENT 0.1 NLO 
calculations (full lines) for four out of eight investigated bins in Q. The spectra for 
(Q) = 7.5-57.9 GeV are multiplied by factors of 10 n , n = 0,1,2,3. The error bars 
represent statistical uncertainties only. Note that the first data bin was left out in both 
normalizations. 
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10.1 The Tube Model 

Measurements of the momentum spectra of secondary particles produced by a primary 
parton of energy E$ = Q/2 approximately obey rather simple scaling laws ||. Charac- 
terizing a particular secondary by E, pt, p\ and m with respect to the direction of the 
original parton, the distribution inp t , p{pt), decreases exponentially almost independently 
of Eq and pi, whereas the distribution in pi rises proportional to Eq. Making use of the 
rapidity variable y R (s. eq. ( |5.3|) ), 



1 E + Pl E+ Pl 

y R :=-hx— =ln =, 10.1 

2 E-pi y/pj + rn? 

with the property dy R = dpi/E, this can be expressed in the form of 

*{y R ,Pt) = B • pfa) . (10.2) 

B is the (constant) number of hadrons per unity in rapidity and $ the distribution function 
of light hadrons in the produced jet corresponding to a tube in (y R ,p t )-speice. Following 
refs. |5fj, |69[| , the energy Ej and longitudinal momentum of a jet of length y Rm!aL in 
rapidity can be estimated by: 

2/ fi, max OO 

Ej = J J E ■ Bp(p t )dptdy R w A sinh?/ Rimax « p ijmax w E , (10.3) 
o o 

2/ ii, max OO 

Pi,j = J J Plj ■ Bp(p t )dp t dy R w A(cosh?/ Rimax - 1) w E - A . (10.4) 



Recall that E = m t cosh y^ and p\ = m t sinh y R with the transverse mass 

m t := \/p| + m2 ~ Pt (10.5) 
for light hadrons. A is proportional to the average transverse momentum: 

OO 

\:=B J p t p(pt)dp t . (10.6) 
o 

The power corrections to mean r = 1— thrust and jet mass now become 



(p)P ow ^. (10 . 8 ) 



With B pa 2 and (p t ) ~ 0.3 GeV p9 |, we can estimate A ~ 0.6 GeV. Testing this simple 



ansatz with event shapes in ep DIS is the topic of the next section. 
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10.2 1/Q p -Fits 



The input for all fits performed here and in section |10.4| consists of the corrected data 



means of table |7.1| and the perturbative coefficients given in tables |9.2| - |9.8| together with 
their statistical uncertainties. Systematic uncertainties are not taken into account; they 
will be analyzed in section |10.5| . 

The standard fit procedure employs the method of least squares (s. e.g. J7D[ ) to de- 
termine the best estimates for the parameters of the supplied model. Technically, this is 
accomplished by using the program MINUIT [71] which additionally provides the x 2 - value 



and, where appropriate, the maximally encountered correlation k. In the tables collecting 
the results, the reduced x 2 - va lues, i.e. x 2 divided by the number of degrees of freedom 
(dof), are entered. For one- and two-parameter fits we have seven respectively six dofs 
yielding a maximal x 2 /dof ~ 2 at a confidence level of 5%. 

The first test carried out is to allow a variation of a s (Mz) in the perturbative expres- 
sion 

(F) pert = c 1>F (x) a s (p R ) + (^c 2 ,f(x) + ^ In ^c 1>F (x)j « s 2 (^) (10-9) 



only, where, for the time being, p R is identified with Q. As can be seen from table 10.1 . 
this works surprisingly well except for Bp. Seemingly, the decrease of {F) with Q can 
be accounted for by the logarithmic Q-dependence of a s as long as a s (M z ) may be 
adapted to the corresponding hadronization correction. Compared to the world average 
a s (Mz) = 0.119 0, we obtain high results for Tp and Bp, very high ones for p F , Tc, Cp 
and lower values for yjj, y^ t . Taking this as an indication for the impact of hadroniza- 
tion, it can be concluded that our event shapes are affected by medium, large and small 
(negative) corrections respectively. 

Sometimes, the ratio of NLO to LO results, the i^-factor, is considered to signal big 
effects of higher orders. However, although evaluating them for Tp and pE in Q-bin one 
to be K Tp = 1.28 and K PE = 1.16, we get ratios of the data means with respect to NLO 
of 1.30 and 2.07 calling this argument into question. 

Writing 

(F) = (F) pcrt + (F) pow (10.10) 
according to 0, we extend our fit to include terms of the form 



and 

( i? ) POW = ^2- ( 10 - 12 ) 



At first, a s (Mz) will be kept fixed to 0.119 yielding estimates for A and p as given in 
table |10.1| . With the exception of Bp, the x 2-varues worsen for both formulae, even 
dramatically in the case of p/Q 2 . Simple l/Q 2 -corrections are therefore ruled out, but 
notably y^ t exhibits the best behaviour of all. Also, pure l/Q-terms do not work very 
well, yet, at least the order of magnitude for A ~ 0.5-1.0 GeV is confirmed for Tp and pe, 
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1 T?\ 


oc s {M Z ) 


X /dot 


A/ loreV 


X 2 /dof 


HI VoreV 


X /dot 


(r P ) 


0.1311 ±0.0005 


0.5 


0.705 ±0.028 


1.4 


7.17 ±0.32 


22.8 


(Bp) 


0.1341 ±0.0005 


24.0 


0.549 ±0.016 


0.7 


5.22 ±0.16 


28.0 


(Pe) 


0.1601 ±0.0006 


0.6 


0.539 ±0.007 


38.4 


5.17 ±0.08 


282.5 


(to) 


0.1520 ±0.0005 


1.5 


0.734 ±0.011 


29.2 


7.12 ±0.13 


212.2 


(Cp) 


0.1481 ±0.0004 


5.4 


2.325 ±0.032 


66.2 


20.82 ±0.36 


326.6 


(Vfj) 


0.1110 ±0.0005 


1.7 


-0.185 ±0.012 


2.0 


-1.74 ±0.13 


10.5 


(Vkt) 


0.1082 ±0.0007 


1.8 


-0.360 ±0.021 


1.4 


-4.21 ±0.27 


9.3 



Table 10.1: Results of one-parameter fits according to eqs. ( J1U.9 ) — ( |10.12 ) for the event 
shape means. Uncertainties are statistical only. 



although a relative factor of two is not found between them. According to our previous 
expectations, Cp receives a larger A, whereas for y$j and y^ t negative values are assumed. 

In summary, this suggests that logarithmic as well as power-like Q-dependences play 
a role. The idea to fit both, a s (M z ) and A, however, does not lead to the desired outcome 



as shown in table |10.2j . Due to the extreme anti-correlation found between them, there is 
a tendency to minimize the power contribution at the cost of a very high a s {Mz)- Merely 
Bp and the y event shapes produce correlated but reasonable numbers. Triggered by the 
x dependence of (F) pcrt , one could assume that (F) pow should be a function of x as well. 
Yet, for corrections with power p = 1 this is not expected [f72| . For an improvement a 
more sophisticated model is needed. 

At last, an explicit check on the power p was performed by fitting 

(F) pow = — (10.13) 

\ / Q p \ ) 

keeping again a s (Mz) fixed. The results for p lie between 0.5 and 1.0, but once more v 
and p are strongly correlated. 



10.3 The Model of Dokshitzer, Webber et al. 

For shorthand notation the model introduced in the following is labelled after Dokshitzer 
and Webber, the authors of ref. H, where the essential formula ( |10.15| ) needed here is 
first presented. Note that besides these initiating authors many others contributed to the 
currently very active field of power corrections in QCD, s. e.g. J73J and refs. therein. 

As demonstrated above, pQCD applied in the form of NLO integration programs is not 
sufficient to describe our data. Ignoring the fact that the calculation of still higher orders 
will become forbidding due to the tremendous amount of Feynman diagrams to compute, it 
may not even be a good idea, since the QCD perturbation series is not necessarily expected 
to converge! One known source of singularities are diagrams with a chain of fermion loops 
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/ Z7\ 


oc s {M Z ) 


A/ UeV 


X /dot 


K/ /O 


ol s \M z ) 


HI Uev 


X /dot 


K/ /O 


(tp) 


0.1318 


-0.044 


0.6 


-99 


0.1314 


-0.23 


0.6 


-89 


(Bp) 


0.1195 


0.532 


0.8 


-92 


0.1279 


3.07 


1.7 


-68 


(Pe) 


0.1652 


-0.089 


0.5 


-99 


0.1614 


-0.27 


0.4 


-87 


(rc) 


0.1573 


-0.153 


1.3 


-99 


0.1535 


-0.57 


1.1 


-86 


(Cp) 


0.1609 


-1.393 


2.0 


-99 


0.1519 


-4.84 


1.0 


-85 


(Vfj) 


0.1142 


-0.076 


1.6 


-97 


0.1119 


-0.28 


1.7 


-82 


(Vkt) 


0.1216 


-0.453 


1.6 


-99 


0.1050 


1.33 


1.6 


-93 



Table 10.2: Results of two-parameter fits according to eqs. ( |10.9j )-( fi0.12| ) for the event 
shape means. 



inserted into a gluon line. These renormalon chains lead to a factorial divergence of the 
coefficients in the series expansion ( 1U9D . At first sight, we are at a loss here. But on the 
contrary, one can even gain more insight into non-perturbative effects of QCD by studying 
these divergences [[l6|, [73| . In a manner of speaking, one approaches the non-perturbative 
regime from the perturbative side. Employing this renormalon approach |73|, one ends up 
with ambiguities of the series expansion of the general form (A/n R ) p modulo logarithms 
\ia q ([i R /A), yet the exact result for an observable has to be single- valued, of course. Hence, 
power-like terms have to be added to regain definiteness. 

Practically, this scheme can be exploited to give the leading corrections in terms of 
powers p and q for an observable, but the absolute normalization is unknown. Introducing 
the notion of a universal infrared- finite effective coupling a e g (fi R ) responsible for soft gluon 
emissions, the coefficients for different observables to the same power p can be related 
at the expense of a new but universal non-perturbative parameter a p _i(/i 7 ). Above an 
infrared matching scale with 

A</i I ^2GeV<Q, (10.14) 
this effective coupling has to coincide with the normal perturbative one a s (fi R ). For event 



shapes in e + e annihilation, this concept was first applied in and extended in [24| : 

" / ' ' Po (\ fh* K 1 

Pi Po P 



^pow _ Of / Pi Y / Pr 



P \Pr 



Pi 



oip-xilit) - a 8 (fi R ) - £M In — + ^ + - ) o? s {Pi 

Z7T \ Hi Po 



Here, hp is an F-dependent but calculable coefficient and 



K = N f 

6 2 9 ; 



10.15) 
10.16) 



accounts for differences between the MS and the CM W renormalization scheme |75|, |7g] : 

K 



10.17) 
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The subtractions proportional to a s and a 2 s should be evaluated for the same number of 
flavours, here Nf = 5, as in the perturbative part eq. (|10.9|) and serve to avoid double 
counting. 

In terms of the effective coupling, can be written as the moment 



a 



p-i 



V f 
— / k p 

Fi J 

o 



l a e e(k)dk . 



(10.18) 



i.e. ao (/•*/) essentially corresponds to an average effective coupling at low scales and is 
deduced from e + e~ data on thrust to be ~ 0.5 for \i I = 2 GeV ||. 

For our event shapes, p — 1 except for yu t where p = 2 is expected [f24], |2l| . In 



that case, other mechanisms can give contributions to the same order 1/Q , preventing a 



straightforward calculation of a 2/fe< 
when performing the transition 



The other coefficients can be gathered from 



ap 



16 

67T 



(10.19) 



to normalize to the dp- value for r in e + e~ physics [0. 

In the form of the dispersive approach [77], this technique was applied to event shape 
distributions [78| and other subjects as well [|7j|, 80]. In the meantime, however, it was 



discovered that some ambiguities unfortunately remain in deriving ap- values. Refining 



the investigations to two- loop level |8T]. |82| , the problem could be resolved and universality 
was reconstituted provided an additional factor, the Milan factor 



M 



7.311 - 0.052iV 



/ 



0o 



(10.20) 



is applied. Recently, it was confirmed to hold also for DIS event shapes |p3| . Due to a 



mismatch in definitions compared to ( |10.15[ ), however, we have to multiply by 2/tt ■ Ai 
1.143 where, in general, Nf = 3 should be used p4| . 
Putting everything together, we get for (F) pow 



(F) 



pow 



ap 



32 M 



3tt 2 



P 



2tt 



ln^ + ^+ 1 



Po P 



(10.21) 

with the coefficients given in table |10.3 . 

Note that there is no logarithm with power q any more. Such a term with q = 1 was 
predicted for Bp only. In the previous publications [^, this could not be supported 



and the theoretical calculations for it have been reexamined [76|, leading to the much 
more complicated factor 



7T 



L B P 



2 \l 3« Si cMw(e 



- 3 / 4 Q) 



+ 7- 77 -0.6137 + 0(1) 
4 16 



(10.22) 
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(1 p 


P 


ULp—\yjJjI — Z Vorfc! V ) 




-\/ 2 /Hnf 

X / ao1 


rt/ /0 


(r P ) 


l 


1 


0.480 ± 0.028 t dM2 


0.1174 ±0.0030 ±S:Soax 


0.5 


-97 
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1/2 • a' Bp 


1 


0.491 ± 0.005 to° f 6 


0.1106 ±0.0012 tof 057 


0.7 


-58 


(Pe) 


1/2 


1 


0.561 ± 0.004 toml 


0.1347 ±0.0015 to° lll 


1.2 


±7 


(re) 


1 


1 


0.475 ± 0.003 _o;o4| 


j~\ -1 /™v /"v A 1 /~\ /~\ /~\ ~i A —1—1 1 111 11(1 

0.1284 ± 0.0014 Io;oo92 


1.3 


±19 


(Cp) 


37r/2 


1 


0.425 ± 0.002 to° Q H 


0.1273 ± 0.0009 1q;oo93 


0.9 


±63 


(Bp) 


1/2 


1 


0.602 ± 0.014 to° ll 


0.1164 ±0.0011 tomlt 


0.7 


-84 


(Vfj) 


1 


1 


0.258 ±0.004 


0.1044 ±0.0017 


1.9 


-61 


(Vfj) 


-1/4* 


1 


0.460 ± 0.079 


0.1183 ±0.0027 


1.5 


±99 


(Vkt) 


It 


1 


0.431 ±0.034 


0.0462 ±0.0129 


2.6 


-99 


(Vkt) 


It 


2 


0.586 ±0.182 


0.1047 ±0.0019 


1.6 


±93 


(Vkt) 


-30* 


2 


0.326 ±0.034 


0.1206 ±0.0040 


0.3 


±99 



Table 10.3: Results of two-parameter fits a la Dokshitzer, Webber et al., eqs. ( |10.9|) 
and ( 10.21| ), for the event shape means. In case of Bp, the outcome without the additional 
factor a' Bp , eq. |10.22j , is given as well. The starred coefficients of the y variables are derived 
from a circular fit procedure described in section |10.4.3| , whereas the coefficients marked 
with t's are guesstimated. Where reasonable, statistical and total systematic uncertainties 
are presented. 



10.4 Fits a la Dokshitzer, Webber et al. 



10.4.1 Two-Parameter Fits 



Two-parameter fits according to eqs. (|10.9|) and (|10.21|) are performed for each event 
shape separately. The results are compiled in table |10.3| , where for comparison the fit of 
(Bp) without the new factor a' Bp is given as well. For a discussion of the quoted total 



systematic uncertainties s. section 10.5 



Omitting at first the y variables, one immediately recognizes reduced correlations as 
compared to the tube model — with the exception of rp — and reasonable x 2 /dof-values. 
The five obtained a (fij = 2 GeV)'s scatter around the expectation of ~ 0.5 by about 20% 
and within errors are compatible with it. The uncertainties are dominated by systematics 
of 5%-10% being at least twice as large as the statistical ones. 

Turning to a s (Mz), one observes a two- fold ambiguity. The event shapes employing z* 
as event axis prefer low couplings, whereas the other three lead to a s (Mz) ~ 0.130. Con- 
sidering the errors of l%-3% statistically and 5%-9% systematically, they are consistent 
with each other. Yet, it does evoke some suspicion. 

Figs. |10.3| and |10.4| show the fit results in comparison with the corrected data as well as 
the LO and NLO predictions of DISENT 0.1 with respect to each corresponding a s (Mz). 
The NLO band represents a variation of the renormalization scale of 1/2 < /z|/ GeV 2 < 2, 
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Figure 10.3: Corrected means (full symbols) of r P , B P , p E and r c as a function of Q. The 
error bars represent statistical uncertainties only. The full line corresponds to a power 
correction fit according to eqs. ( |10.9| ) and ( |10.21[ ). For comparison the LO (dashed) and 
NLO (band) predictions of DISENT 0.1 are shown. The band reflects a variation of the 
renormalization scale of 1/2 < /j, 2 R / GeV 2 < 2. 
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Figure 10.4: Corrected means (full symbols) of Cp, yfj and y^ t as a function of Q. The 
error bars represent statistical uncertainties only. The full line corresponds to a power 
correction fit according to eqs. (|10.9|) and ( |10.21| ) except for y^ t where no power contri- 
bution was assumed. For comparison the LO (dashed) and, if appropriate, NLO (band) 
predictions of DISENT 0.1 are shown. The band reflects a variation of the renormalization 
scale of 1/2 < fij GeV 2 < 2. 
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Figure 10.5: Results of the two-parameter fits a la Dokshitzer, Webber et al. for the means 
of Tp, Bp, pe, tc and Cp in the form of la- and 2a-contours in the (a s , ao)-plane for 
statistical uncertainties only. 



s. section |10.5| for details. Each fit taken separately, the data are nicely described. Putting 
all five in terms of (a s , ao)-pairs into the one fig. |10.5| , the discussed discrepancies become 
obvious. But recall again that systematic uncertainties are treated individually and there- 
fore only statistical uncertainties are included! 

With yfj, one obtains a good fit but rather low numbers for ao(/-t/ = 2 GeV) and 
a s (Mz), leading to the conclusion that the coefficient a VfJ given in the proceedings ]24| 
may not be appropriate. Referring back to section 10.2 , yjj does exhibit smaller hadroni- 
zation corrections than thrust. The curve in fig. |10.4| was produced with a yfJ = — 1/4 
derived in section |10.4.cl . 

Concerning y^ t , the coefficient a Vkt is not known and the 1^ in table |10.3| is pure 
guesswork. The fit for p = 1 works not very well supporting the conjecture of a l/Q 2 - 
term instead of 1/Q. In fact, both results obtained with p = 2 look more reasonable. 
Nevertheless, fig. |10.4| contains the outcome according to eq. (|10.9| ) without any power 
correction. 



10.4.2 Combined Fits 

Based on QCD, the universality of a s is on firm theoretical grounds. Not so, however, 
ao since here additional assumptions are made. It was therefore checked if fits allowing 
for different «o's are possible while enforcing one a s [Mz)- Correlations between the 
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OOII1D. ncs. 


TP) -DP 


Pe, t c Cp 


TP; PP; r C; ^P 


Tp; -Dp, pp, Tc, Op 


tto(Tp) 


0.532 ±0.012 


— 


0.374 ±0.010 


0.414 ±0.009 


a (B P ) 


0.488 ±0.005 






0.443 ± 0.006 






0.561 ±0.004 


0.561 ±0.004 


0.562 ±0.004 


ao(Tc) 




0.476 ± 0.003 


0.476 ±0.003 


0.474 ± 0.003 


a (Cp) 




0.428 ±0.002 


0.427 ±0.002 


0.421 ±0.002 


« S (M Z ) 


0.1115 ±0.0011 


0.1291 ±0.0007 


0.1285 ±0.0007 


0.1245 ±0.0006 


xVdof 


11.6/13 


38.4/20 


55.5/27 


222/34 



Table 10.4: Results of combined fits a la Dokshitzer, Webber et al. allowing for a separate 
ctoiPi — 2 GeV) for each event shape, but enforcing one a s (Mz). Several combinations 
are tried; uncertainties are statistical only. 



observables were ignored for that purpose. According to the systematic behaviour found 
previously concerning the a s (Mz)-results, several combinations are tried. The outcome 
is collected in table 110.4. 



The first two arrangements reflect the separation in Tp, Bp on the one hand, and 
Pe, tc, Cp on the other. As expected, the numbers show resemblance to what we have 
obtained before. Yet, including step by step first r P and then B P into the p#, r c , Cp 
triple, the x 2 /dof -values worsen considerably.^ Both times the fit results are dominated 
by the latter, which in addition to be more numerous, also have the smaller statistical 
uncertainties. a s [Mz) is forced to remain high, and thereby a (/i/ = 2 GeV) for t p and 
Bp is shifted. 

One concludes that combined fits with one a s (Mz) are possible, but notwithstanding 
the freedom to choose an individual «o for each event shape, they suffer from systematic 
effects. In order to determine a s (Mz), it would be interesting to know what their origin 
is and what biases are implied. 



10.4.3 a F -Fits 

Motivated by the problems encountered with the a^-coefficients for the y event shapes, 
one can try to do a three-parameter fit of ao, a s (Mz) and dp. However, this set of 
parameters is strongly correlated and with eight data points only it does not properly 
converge for Tp, Bp, pp, Tc or Cp. The y variables are somewhat of an exception here, 
giving a (pj = 2 GeV) = 0.58 ± 0.20, a s (M z ) = 0.134 ± 0.015 and a VfJ = -0.71 ± 0.18 
at x 2 /dof = 1.6 for y fJ , and a^p, = 2 GeV) = 0.35 ± 0.05, a s (M z ) = 0.123 ± 0.006 and 
a Vk t = —33.2 ±3.0 at x 2 /dof = 0.23 for y kt respectively. Despite looking quite reasonable, 
one nevertheless has to exercise caution because of the very strong correlations. 

1 Note that with increasing dots from 6 for the two-parameter fits to 34 in the last combination, the 
limit for the 5% confidence level falls below 1.5. 
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Table 10.5: Results for fits of the coefficients while presupposing the validity of eq. |10.21 
with Q!o(a*i = 2 GeV) = 0.5 and a s (Mz) = 0.119. Uncertainties are statistical only. 



Being constrained that way to less than three parameters, we start from the other side 



by deriving ap while presupposing the validity of eq. |10.21| with cko(a*j = 2 GeV) = 0.5 and 
a s (Mz) = 0.119. Table |10.5| compiles the determined coefficients, which, although some 
rather large ^ /dof's occur, are of a similar magnitude as the theoretical ones. Notably, 
Ufj gives ap ps —1/4 instead of 1! 

As a consistency check, this new value was reinserted to produce the two-parameter 
fit in fig. [T01] with 5 (/i/ = 2 GeV) = 0.46 and a s (M z ) = 0.118, s. again table (TTT^. 



Performing the same, i.e. p = 1, circular procedure with y kt , one gets ao(fJ-i = 2 GeV) = 
0.34 and a s [Mz) = 0.112 which looks much better than the entry in table |10.3| , but still 
is not satisfactory. Since ai is basically unknown, we must determine it in addition to 
ap in case of p — 2. Refitting then with ap = —30, one obtains ai(/i 7 = 2 GeV) = 0.33 
and a s {Mz) = 0.121, s. table |10.3| . Both refits, however, exhibited strong correlations 



between a respectively ai and a s (M z ). 

In a last test, all coefficients were multiplied by factors of 2 respectively 1/2 to study 
the effect on the a^/ 1 / — 2 GeV), a s (M z ) results. Having fig. |10.5| in mind, the outcome 
corresponds approximately to shifts along the main diagonal with large ap's inducing 
smaller o7 (/U 7 = 2GeV)'s and a s (M 2 )'s. 



10.5 Evaluation of Systematic Uncertainties 



The general procedure followed here is to repeat the fits for a variation in every prominent 
origin of systematic effects. The obtained discrepancy compared to the standard result 
is attributed to a corresponding uncertainty. In case of two deviations in one direction 
for the same primary source, e.g. an up- and downwards modification of an energy scale, 
only the larger one is considered for the evaluation of the total uncertainty. The latter 
is derived from all contributions by adding them up quadratically. An exception is the 
unfolding, whose influence is estimated in the same way as explained in section [775] on the 
final results of the data means. 
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Altogether, the following studies have been performed to estimate systematic effects: 

• Experimental uncertainties: 

1. Usage of four unfolding procedures 

2. Variation of the electromagnetic energy scale of the calorimeters by ±1%, 2%, 
3% depending on z irnp 

3. Variation of the hadronic energy scale of the LAr calorimeter by ±4% 

• Theoretical uncertainties: 

1. Variation of the renormalization scale /i| = Q 2 by factors of 2 and 1/2 

2. Variation of the factorization scale /i| = Q 2 by factors of 4 and 1/4 

3. Variation of the infrared matching scale p l by ±1/2 GeV 

4. Selection of the MRSA'-105 and MRSA'-130 pdfs @ with lower respectively 
higher intrinsic a s (Mz) than in the standard pdfs 

5. Usage of the completely different pdf set CTEQ4A2 |35| with similar a s (Mz) 
The experimental sources have already been described in section [7.5|, but concerning 



the theoretical ones, some remarks are in order. The renormalization scale fx R as well as 



the factorization scale /i F (s. [50] for details) are arbitrary in the sense that in a complete 
theory the calculations are not allowed to vary with any specific choice. Yet, in reality we 
have only an approximative theory at our disposal yielding residual dependences due to the 
neglected higher orders. To avoid the appearance of large logarithms in the computations, 
it is recommended to identify them with a process relevant scale which for our purposes 
is always Q. To estimate the effect of omitted higher orders, it is conventional to vary 
them by an arbitrary factor of 4 concerning the squares fj, 2 R , fj, 2 F . In the case of fi 2 R , we had 



to reduce this factor to 2 because of the condition (|1(J.14[): A <C fij m 2 GeV <C Q. Taking 



this into account, varied results can be achieved quickly via eq. (|8.21|) . With respect to 
fi F , however, the complete calculations have to be redone. 

The change of the infrared matching scale /i 7 by ±1/2 follows from the values used in 
the original proposal |J. Note that this contributes only to the uncertainty for a s (M z )\ 
ao explicitly depends on fi r . 

The next point accounts for the fact that implicitly a s (Mz) has already been used in 
deriving the pdfs from data which in turn may bias our computations. The same is true 
for the choice of parameterization of the pdfs. For those reasons, three alternative sets, 
two with different assumptions on a s [Mz) and one with approximately the same a s (Mz) 
but another parameterization, have been selected for a reevaluation with DISENT 0.1. 

All uncertainties described are presented graphically in figs. |10.6| - |T0^8i separately for 



each of the five event shapes where the coefficients ap are known. Without such a pre- 
diction for yfj and yk t , a study of systematic effects is omitted. 

Inspecting the figs. |10.6| - |10.8l again, it stands out that Tp and B P behave differently 



compared to pp, t c and Cp. For the latter, the systematic uncertainties are clearly 
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Figure 10.6: Systematic uncertainties of ao(/J>i = 2 GeV) (left) and a s [Mz) (right) for rp 
(top) and Bp (bottom). 
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Figure 10.7: Systematic uncertainties of ao(/U/ = 2 GeV) (left) and a s [Mz) (right) for pe 
(top) and Tc (bottom). 
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dominated by the variation of the renormalization scale. In case of rp and Bp, this is 
much less pronounced. The larger influence of experimental uncertainties probably can be 
related to the explicit reference to the boson axis implied in the definitions of Tp and Bp. 



10.6 Cross-Checks 

Concerning systematic effects some additional cross-checks were performed: 

1. In order to ensure a reliable measurement of the scattered electron, partially in- 
efficient regions between calorimeter modules (0-cracks) or wheels (^-cracks) are 
avoided. This selectively diminishes the contribution to certain phase space regions 
and is corrected for by the procedure described in section [7]2|. Nevertheless, it was 
checked that the actual influence is negligible, even without unfolding. 

2. The derivation of the power correction coefficients does not account for phase space 
constraints such as the cuts nos. ^| and § imposed on the scattered electron. Expe- 
rimentally, they are absolutely necessary, but for testing purposes the NLO calcula- 
tions were repeated without these cuts. Except for (yk t ), which increased by about 
8% in the low Q 2 region, all other mean values changed by less than 2%. 



Chapter 11 

Summary and Outlook 



The first measurement of event shapes in ep DIS [23|], covering a range in Q from 7 GeV 
up to 100 GeV in a single experiment, has been substantially improved and extended to 
include the additional variables p#, Cp, yfj and as well as new data taken in 1997. The 
available statistics in the high Q 2 sample was more than doubled, facilitating a bisection 
of the highest Q-bin. The improvements encompass i.a. a more precise calibration of the 
energy of the scattered electron, an unfolding of the data distributions applying a Bayesian 
approach and a full account of radiative QED corrections. Taking into consideration the 
predominant sources of systematic influences, i.e. the accuracy of the electromagnetic 
and hadronic energy scale of the calorimeters and the unfolding procedure, corrected 
mean values of the event shape spectra together with estimates on their statistical and 
systematic uncertainty could be derived. 

All event shape means exhibit a strong Q-dependence. They decrease with rising Q, 
i.e. the energy flow becomes more collimated. For the event shapes defined in the current 
hemisphere alone, it is interesting to compare with corresponding ones in e + e~ collisions. 
This was already done in for (tq) and (rp) as well as (pq). Despite the fact that the 
QCD dynamics in ep DIS is differentia large similarity not only in the Q-dependence, but 
also in absolute value to their e + e~ counterparts could be established. When consulting 



reanalyzed data of the JADE Collaboration [87] , this still holds for the C parameter, but 



not for the jet broadening. Bp is about a factor two larger than the wide as well as the 
total jet broadening. 

Employing calculations in pQCD to NLO accuracy, the influence of parton fragmenta- 
tion and hadronization could be studied. As an improvement, the approximative treat- 
ment of the x-dependence of the perturbative coefficients applied in |23j has been replaced 
by a more precise bin-wise evaluation. A direct comparison with data revealed that the 
investigated event shapes can roughly be categorized into three classes: 



1 There are neither BGF diagrams in e + e physics, nor is there a predefined event axis. In addition, 
it is impossible to have an empty hemisphere. 
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1. The event shapes Tp and Bp employing the z*-axis as event axis are sensitive to 
QED radiation. They show hadronization corrections of medium size. 

2- Pei t c an d Cp are strongly affected by the non-perturbative transition to hadrons. 
Their by far largest systematic uncertainty is due to the variation of the renormali- 
zation scale p R . 

3. Both event shapes based on jet algorithms, yfj and yk t , exhibit small (even negative) 
hadronization corrections. 

Fixing a s (Mz) to 0.119, the current world average 0], simple power corrections X/Q 
and p/Q 2 to parameterize the discrepancies between data and theory were tested. For 
A, reasonable results of 0.7 GeV and 0.5 GeV could be obtained for rp and Bp, whereas 
yfj and yk t yielded negative values of —0.2 GeV and —0.4 GeV respectively. For the other 
event shapes the fit delivered poor x 2 /dof's. This is generally the case for p/Q 2 '-terms; 
except for maybe yfj and y^ t , they are ruled out. Fitting a s (M z ) in addition to A or p 
did not help either due to very large correlations. 

Invoking the approach for power corrections initiated by Dokshitzer and Webber, 
the situation improves significantly. For t p , B p , p E , tq and Cp, where predictions for 
the coefficients ap are available, correlations were reduced - - although not as far as 
desirable — and the fits converged nicely. Note that for Bp the new coefficient a' B 
replacing the formerly predicted logarithmic term had to be included. Within 20%, all 
determined ao's are compatible with a universal value of ao{pz — 2 GeV) ~ 0.5. Very 
similar results are achieved with this model when applied to e + e~ data |87], except 
for the wide jet broadening. 

Concerning a s (Mz), an inconvenient spread ranging from 0.111-0.135 is observed. 
Although all five values turned out to be consistent with each other, a simple averaging 
seems to be inappropriate since the outcomes for a s {Mz) again reflect two of the categories 
defined above. Tp and Bp mark the lower bound of the range in a s [Mz) as opposed to 
Pe, tc an d Cp at the upper end. Combining several of these event shapes, the fits did 
work, but their quality severely declined, even allowing an individual a for each.0 As 
long as this systematic behaviour is not better understood, a precise determination of 
a s (Mz) is not feasible that way. Derived from the spread above, it can be stated at best 
that a s (M z ) = 0.123 ±0.012. 

Returning to the y event shapes, the conjectured coefficient of a VfJ = 1 has been 
excluded. Derived from data keeping ao(/-t/ = 2 GeV) = 0.5 and a s (Mz) = 0.119 fixed, 
a value of about —1/4 can be predicted. Without any coefficient given for y^ t , it can 
only be checked for the power p of the 1/Q P term. Due to the smallness of the observed 
hadronization corrections in comparison with experimental uncertainties, however, it is 
very difficult to extract them from data. In fact, the performed investigations hint at 
powers p > 1, but they are not conclusive enough to claim more. 



2 Note that correlations were ignored. 



117 



The achieved measurements can be complimented in several ways: 

• Instead of looking for the Q-dependence alone, the data could additionally be divided 
into bins in x. This is done by the ZEUS Collaboration [^6] . 

• In e + e~ annihilation the model of Dokshitzer and Webber et al. has been extended 



to differential distributions [78, 76]. Although no predictions for ep DIS exist, it 
is straightforward to just look what comes out, since the unfolded event shape 
spectra are at hand. Actually, a first check was already performed f49| and revealed 
inconsistencies even between fits of the mean values and distributions of one shape 
variable. Resummed calculations which are not available for ep DIS may be helpful. 

By further subdividing a distribution into intervals as done by the DELPHI Collabo- 



ration, it is possible to examine where specific power-like contributions occur 3S, 90 



Alternatively, higher moments of already investigated event shapes can be studied. 



First steps into that direction were taken by the OPAL Collaboration [BT 



• Another obvious extension is the usage of new variables, possibly with small ha- 
dronization corrections. 

Especially concerning the last point, however, the question arises of what to do with 
them. For y^ t it was already noticed that with the given experimental uncertainties power 
contributions can not be reliably extracted. Yet, just in this case it could be interesting to 
perform direct fits of a s [Mz) with the distributions at high Q 2 . Remaining hadronization 
corrections must then be estimated by the usual procedure of correcting to the parton 
level of MC models. At the same time, this could be done for the event shapes of this 
analysis, such that a comparison of the fit results in a s (Mz) could provide insight into 
systematic effects in one (or both) of these methods. 

At last, it should be mentioned that the per se sharp angular cut-off of 90° in the Breit 



frame is neither theoretically nor experimentally unproblematic, (s. sections [4.2| and |5.5| ). 
It has been suggested to replace it by a kind of weighting procedure to ensure a smooth 
transition. On the other hand, one could try to separate the remnant by employing e.g. 
the fct-algorithm and treat then all remaining objects as input for the derivation of event 
shapes. 
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